GEOMETRY OF THE PFAFF LATTICES 
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Abstract. The (semi-infinite) Pfaff lattice was introduced by Adler and van Moerbeke [2] to 
describe the partition functions for the random matrix models of GOE and GSE type. The 
partition functions of those matrix models are given by the Pfaffians of certain skew-symmetric 
matrices called the moment matrices, and they are the r-functions of the Pfaff lattice. In this paper, 
we study a finite version of the Pfaff lattice equation as a Hamiltonian system. In particular, we 
prove the complete integrability in the sense of Arnold-Liouville, and using a moment map, we 
describe the real isospectral varieties of the Pfaff lattice. The image of the moment map is a convex 
polytope whose vertices are identified as the fixed points of the flow generated by the Pfaff lattice. 
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1. Introduction 

This paper centers around the finite version of the Pfaff lattice hierarchy. The Pfaff lattice is 
a Hamiltonian system, and has many connections to lattices of Toda type. Our goal is to show 
the complete integrability of the Pfaff lattice, and to describe the geometry of the iso-level sets 
determined by the integrals of the lattice. 



1.1. Lie algebra splitting related to SR-factorization. We introduce the Pfaff lattice by start- 
ing with a Lie algebra splitting (see [TJ [5] for the semi- infinite case) , 

sl(2n,R) = £®sp(n,R) . 

Here the splitting is defined by the projections of an element X £ sl(2n, R), 

ixt{X) = X- — J (X+) T J + ~(X - JXfj) 

tt sp (X) =X+ + J (X+f J+l(X + JX*J) 

where Xq is the 2x2 block diagonal part of X, X + (resp. X_) is the 2x2 block upper (resp. lower) 
triangular part of X, and J is the 2n x 2n matrix with the 2x2 diagonal blocks, 




We denote sp(n) = {X £ sl(2n) : a(X) = X} with the involution a(X) := JX T J. One finds that 
dimsp(n) = 2n 2 + n. Each element Y £ t is expressed by a 2 x 2 block lower triangular matrix with 
2x2 block diagonal part given by Yq =diag2(fci/2, ■ • • , knh) with Y^j=i kj = 0' an d h represents 
the 2x2 identity matrix. We also use the 2x2 block decomposition for g = sl(2n), 

= 9+ © flo ffi 0- = g>o © 0- = 0+ © s<o , 

where g± represents the set of 2 x 2 block upper (lower) triangular matrices, and go the set of 
2x2 block diagonal matrices. With the inner product on s[(2n) defined by (X,Y) := tr(XY), one 
can define the dual spaces 6* = sp(n) ± and sp(n)* = 6 . With the identification g = g*, we have 
s\(2n) = B* © sp(n)* with the projections, 

Tr e , (X) =X+-J (X+f J + i(X - JXq J) 

vr sp , (X) = X- + J (X+f J + ~(X + JX^J) . 

Notice that both t,?- 1 - C 0<o ■ Also t* = sp(n) ± = {X £ s((2n) : <r(X) — -X}, which gives a vector 
space decomposition based on the involution a(X) = JX T J, that is, sl(2n) = sp(n) 0sp(n) _L . 
The sp(n)* is the set of 2 x 2 block lower triangular matrices with block diagonal part of the form 
diag(ai, — ai, a2, —a-2, ■ ■ ■ , a n , —a n )- Note that dim(sp(n)*) = dim(sp(rt)) = 2n 2 + n. 

Remark 1.1. The splitting of Lie algebras used here is induced by the SR-factorization of SL(2n, M) 
which is an Sp(n, R) version of QR-factorization. This factorization states that for a dense set of 
A £ SL(2n,R) we can find Q £ G e with Lie(G t ) = t and P £ Sp(n,R) such that A = Q~ 1 P (see 
Theorem 3.8 in [4]). However Theorem 3.7 in [4] shows that the set of A £ SL(2n,C) such that the 
SR-factorization exists is not dense, for this reason we restrict our analysis to real matrices. Remark 
3.9 in [4] states that the SR-decomposition as defined in that article is not unique; however our 
choice of Gt gives a unique decomposition. 
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Let us now define a matrix L £ s[(2n) , 

L = £ + k with £ e Sp(n)* , k G fi* . 
In particular, we choose k to be k = diag 2 (/3i/2, ■ • ■ , Pnh) + e with X)J=i Pj = an d 





/ o 2 
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Here I 2 is the 2x2 identity matrix, 2 is the 2x2 zero matrix and e_ = ^ J . Then the matrix 

L is a lower Hessenberg matrix, and the number of variables in L is given by 2n 2 + 3n — 2. For 
example, L g sl(4) can be written in the form, 



(1.1) 



L = 
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The number of variables is 2 - 2 2 + 3- 2- 2 = 12. 

This choice of the matrix L corresponds to the minimal nontrivial example of the Pfaff lattice 
considered in this paper. 

1.2. Hamiltonian structure. We define a Lie-Poisson bracket {F, H} S *(L) for functions F, H on 
g* using the i?-matrix R = \{tti — 7r sp ) with the splitting g = t®sp(n): For any I,Yeg, we define 
the i?-bracket 

[X,Y] g = [RX,Y] + [X,RY] 

= [ir t X,n t Y}-[ir sp X,Tr sp Y}. 
Then the Lie-Poisson bracket is given by 

(1.2) {F,H} g ,(L) = (L,[VF,VH] g ). 

The Pfaff lattice for L e g* is defined by a Hamilton equation with Hamiltonian function H , i.e. 



f t ={H,L} A L) 



(1.3) 

(1.4) = ir sr [ir sp VH, L] - 7T t . [TreVff, L] . 

The solution of the Pfaff lattice can be expressed as a coadjoint orbit Ad*^ through L° = £° + n° 

with £° e sp(n)* , k° e and g(i) = exp(iVi/). In general, the coadjoint orbit through L° can be 
written in the form, 

(1.5) L(t) = Ad* g{t) (L°) = t: sp , (PeP- 1 ) + Trt. (Q^Q- 1 ) , 

where g(i) = Q~ 1 P with P e Sp(7i,IR) and Q G Gf. Thus the coadjoint action factors into an 
Sp-action on sp(n)* and a G{-action on £*. 

If the Hamiltonian H is Sp-invariant (i.e. (L, [V-ff, X]) = 0, VX e sp(n)), one can write dL/dt = 
— [7T{Vi/, L]. With the Hamiltonian = -£-Lj-tr(L fe+1 ) (the Chevalley invariants), we define the 
Pfaff lattice hierarchy, 



(1.6) 



8T 

— =-[B k ,L] with B k = TTt(L k ) . 

Otk 



k = 1, . . . ,2n- 1. 
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Then the solutions to ()1.6|) with initial condition L(Q) evolve on the adjoint orbit of both Gt and 
Sp(n,M) through L(Q), that is, ()1.5j) becomes 

L(t) = Ad Q{t) L(Q) = Ad P{t) L(0) with e e ^ L ^ = Q{t)- x P{t) , 

where 9(t, L(Q)) = Ym^L i tkL(0) k (we will discuss more details in Section[3]). 

We now give a canonical Lax matrix L with a form of the lower Hessenberg matrix (jl.ip in the 
sense that we normalize the invariant elements in L under the Pfaff flow of (|1.6|) . Those invariants 
are given by the following: 



Lemma 1.1. For B e t and L a lower Hessenberg matrix, we have 

2i-l 



Proof. We first have 

[-Sj-^]a»-l,2i = (-SIf)2i-l,2i - (£-B)2i-l,2i 



2i-l,2i 



2i-l,2i 



-£'2i-l,2jB2i J 2i , 



where we have used the upper triangular shape of B and the lower Hessenberg form of L. One then 
notes that Bzi-x^i-i = Bii,n from the structure of t and the result follows. The proof of the other 
equations is a similar direct calculation which is left to the reader. □ 

Lemma 1 1 . 1 1 implies that £21-1, 2« > for i = 1, 2, . . . , n, and the traces of the odd principal parts of 
L, Y^j=i Ljj , for % = 1, 2, ... n, are invariants of the Pfaff flow ()1.6j) . We choose values for these 
expressions so that the phase space of the Pfaff lattice is 



L G 0[*(2n,M) : L is lower Hessenberg, L2<-i, 



2i-l 

1, and ^ Ljj 

i=i 



The space with the Lie-Poisson bracket (|1.2|) defines a Poisson manifold, and the dimension is 
given by 

dim(Z R ) = 4 n ( n ~ X ) + n + 2(„ - 1) = 2n 2 + n - 2 . 

The form of L is a finite version of the Pfaff lattice considered in [2] (see (0.12) in Theorem 0.1). 
For example in the case of n = 3, we have 



(1.7) 
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It should be noted that L of this form is given by a point on a Gt-orbit through the companion 
matrix C 7 , that is, L — 'i'C 1 x S~ 1 for some VP G Gt and 



/ 




(1.8) 
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where the 7fc's are the symmetric functions of the eigenvalues, i.e. the characteristic polynomial 
4>2n(z) '■= det(i — zl) = z 2n + X)fc™2( — l) fc 7fe^ 2rl ~ fc - We also denote the real isospectral variety of 
the Pfaff lattice with 7^ £ R as 

2n 

(1.9) Z»( 7 ):={LeZ E :f r , (i)=7r for r = 2, 3, . . . , 2n} = f| {F^{ lr )} , 

r=1 

where F r ^(L) are the invariant polynomials defined by det(L — XI) — z 2n + ^^" 2 (— l) r i 7 ' r .o(L)2 2 ™ _r 
(see (|2?2|)). Note that the set {F rfl (L) : r — 2, . . . . 2n} are equivalent to {H k {L) : k = 1, . . . , 2n- 1}. 
The main purpose of this paper is to discuss the geometric structure of Z^(j) including a complete 
integrability of the Pfaff lattice. 

1.3. Outline of the paper. The paper is organized as follow: 

In Section [21 we discuss the integrability of the Pfaff lattice. While the Chevalley invariants 

H k = -^tr (L k+1 ) , fc = l,2,...,2n-l, 

are Hamiltonians of the lattice, they do not provide enough integrals. We show that the additional 
integrals are given by the generating function 

F L (x, y, z) = dct[(x — y)L + yJL T J — zl] . 

Note that F^(1,0, z) = det(L — zl) = Yi^=i( z j — z )> which generates integrals equivalent to the 
Hk- In addition we note that Fl is Sp- invariant, i.e. Fp^p-i = Fl for P e Sp(n). We follow the 
method of _8j to construct the Casimirs and action-angle variables, with action given by the integrals 
generated by Fi,(x,y, z). We then define an extended Jacobian associated with the spectral curve 
C(L) = {Fl(x, y, z) — 0} which has the genus g(C) = 2n 2 — 4n + 1. The extended Jacobian is given 
by the Abel-Jacobi map of a vector for which 17(C) entries are holomorphic one- forms on C(L) and 
2n — 1 entries are meromorphic one- forms on C(L) (it is the presence of meromorphic one forms that 
makes this Jacobian extended). 

This integrability result is parallel to that of [5] in the sense that the spectral curve Fl(x, y, z) = 
gives the actions (and Casimirs), while an extended Jacobian on the spectral curve gives the angles. 
A fundamental difference in the case of the Pfaff lattice is that the spectral curve F^x^y, z) = is 
singular (Proposition 12 . 1| . The singularities provide n functional relations among the coefficients in 
FL(x,y, z). We then prove that the Pfaff lattice is a complete integrable Hamiltonian system, with 
n 2 — 1 Hamiltonians and n Casimirs, in the sense of Liouville- Arnold fThcorcm l2.4p . 

In Section [31 we review the construction of solutions to the Pfaff lattice by SR-factorization of 
e tL(o)^ rp ne SR_f ac t or i za tion decomposes e tL ^ = Q~ X P where Q £ Gj and P £ Sp(n). A dense set of 
matrices in SL(2n, R) possess such a factorization [4,. Solutions to the Pfaff lattice are then given by 
L(t) = PL(Q)P^ 1 = QL(0)Q _1 . The SR-factorization leads to a Cholesky-type factorization (skew 
Borel decomposition) of the skew symmetric moment matrix M — e tL ^ Je tL ^ — Q^ 1 JQ~ T . With 
the embedding L — * M the Pfaff flow is linear on the space of skew symmetric matrices. We then 
define the T2fc-function as the Pfaffian of the 2fcx2fc principal part of M, and the set (T2, T4, . . . , T2„_2) 
generates the solutions of the Pfaff lattice. We will conclude this section by working out expressions 
for the moment matrix and r-functions which correspond to discrete measure versions of the GOE 
and GSE random matrix models. While the GOE model is a generic case of the Pfaff lattice, the 
GSE model corresponds to a degenerate case of the Pfaff lattice with doubly degenerate eigenvalues. 

In Section [4] we consider generic solutions of the real Pfaff lattice. In particular, we compute 
the fixed points of the Pfaff lattice flow. These are given in the standard way by computing the 
asymptotic limits as t± — > ±00. The result is a sorting property similar to that found for the 
asymptotic limits of the Toda lattice hierarchy. In our case the matrix converges to a 2 x 2 block 
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diagonal structure which sorts the eigenvalues by pairs from smallest pair to largest pair as t\ — > — oo 
and largest to smallest for t\ — > oo (Theorem 14. We then derive the explicit form of the fixed 
points of the Pfaff lattice (Theorem l4.2[) . The total number of fixed points is shown to be (2n)!/2™, 
which is the order of the quotient S% n jWp y where Wp is the parabolic subgroup of the symmetric 
group S2n associated with P2x2- In the case of the GSE-Pfaff lattice, the number of fixed points is 
given by n!, the order of S n . 

We conclude this paper, in Section O by constructing the moment polytope to characterizes the 
isospectral variety of the Pfaff lattice. We first note that the T2k can be considered as a point in 
P(A 2fc K 2 ™). We then define the moment map /Lt : P(A 2fc M 2 ™) — > f)jjj and consider the moment polytope 
as the image of the moment map to describe in detail the geometry of the isospectral variety. We 
show that in the generic case, the image of the moment map is a convex polytope described by the 
irreducible representation of SL(2n,R) on the product P(A 2 M 2 ") x P(A 4 K 2 ") x • • • x P(A 2n ~ 2 R 2 ™). 
The vertices of the polytope correspond to the fixed points of the Pfaff lattice flow, and an open dense 
subset of the polytope is given by the generic orbits of the flow fTheorem l5.ip . The boundaries of the 
polytope consist of the nongeneric orbits of the Pfaff lattice which include the solutions of the KP 
equation. In the case of the GSE-Pfaff lattice, the moment polytope is given by the permutahedron 
of the symmetric group S n , that is, the weight polytope associated to the full flag manifold of 
SL(n, K) (Theorem I5.3j) . This situation is similar to the case of the Toda lattice (see [6], [21"]). 

2. Integrability 

In this section, we prove that the Pfaff lattice hierarchy (|1.6p is a completely integrable Hamilton- 
ian system. Deift, Li, and Tomei [5] consider the case of the generalized Toda lattice, given by the 
Lie algebra splitting sl(n) — b_ o(n), where b_ is the Borel subalgebra of lower triangular matrices 
with zero trace, and o(n) is the set of skew symmetric matrices. They showed the Liouville- Arnold 
integrability of the Toda lattice by finding the full family of Casimirs and action integrals, and the 
conjugate angles. While some of the integrals were found by considering the chops (as in the case 
of the full Kostant-Toda lattice [9]), Deift, Li, and Tomei also introduced the spectral curve 

Fl ° da (x, y, z) = det [(a; - y)L + yL T - zl] = 0, 

for the Lax matrix L G sl(n)* , which is invariant under the S'0(n)-adjoint action, i.e. Fj£^T — F r [ oda 
with q G SO{n). The coefficients of this homogeneous polynomial give the additional Hamiltonians 
and Casimirs needed to find integrability. In addition, the angles conjugate to the actions are a 
section of an extended Jacobian of this spectral curve. 

Our proof of integrability for the Pfaff lattice equations closely follows that in [8] , but with several 
unique structures. The main complications arise from the presence of n nodal singular points on the 
spectral curve. 

2.1. The integrals F r ^(L). Let us consider the following function, 
(2.1) F L (x, y, z) := det[(x - y)L + yJL T J - zl] , 

which defines an algebraic curve on CP 2 with the homogeneous coordinates (x : y : z). We denote 
the curve as 

C(L) := {{x:y:z)e CP 2 : F L (x,y, z) = 0} . 
Note that the function Fl is Sp- invariant, 



Fplp-1 =F L VP G Sp(n) 



GEOMETRY OF THE PFAFF LATTICES 



7 



which also implies (L, [VFl,X]) = for any X E sp(n). Since Fl is invariant under the involution 
i : (x, y, z) —> (x,x — y, —z), Fl can be expanded in the form, 

2n [r/2] 

(2-2) F L (x, y,z) = J2Y, Fr,k(LW (rik) (x, y, z) . 

Here the functions ^ r,k ^ (x, y, z) are homogeneous polynomials of (x, y, z) of degree 2n and are 
invariant under the involution. They are defined by 

<^ 2r > fc ) {x, y, z) = x^-V (y(x - y)) k z 2 ^- r \ 
^ 2r+l - k \x,y,z) = x 2{r - k \2y~x){y{x-y)) k z 2{ - n -^- 1 . 

Then because F is Sp-invariant, the coefficients F r ^{L) provide the integrals of the Pfaff lattice. 
The total number of F r> k{L) is n 2 + 2n — 1 (note Fo,o(L) = 1 and F\fl{L) = are not counted). 
In particular, the functions F Ti q arc equivalent to the Chevalley invariants which are symmetric 
polynomials of the eigenvalues of L, i.e. Fl(1, 0, z) — det(£ — zl). Other invariants F r> k with k ^ 
provide a foliation of the isospectral variety Zyfcy) of ()1.9|) . which will be further discussed in Section 
13.21 To determine a complete picture of the foliation, one needs to find all the Casimirs which define 
the symplectic leaves (symplectic foliation of the Poisson manifold Z^(^)). 

In order to find the independent set of the integrals, we first note that the curve (|2.ip is singular 
at (0,y, z): 

Proposition 2.1. The function Fi,(0,y, z) is reducible, and is given by 

F L (0, y, z) = [pi(yJL + yL T J + zJ)\ 2 , 

where pf (M) indicates the Pfaffian of a skew- symmetric matrix M . 

Proof. First note that the matrix in the Pfaffian, 

M := yJL + yL T J + zJ = -J(-yL + yJL T J - zl) , 

is skew-symmetric. The determinant of a skew symmetric 2n x 2n matrix is the square of the Pfaffian, 
this then proves the Proposition. □ 

Then from Lemma |2.1| we have n functional relations among F r [ r /2](L): 

Corollary 2.1. Each function i 7 '„+fc ! [( n +fc) /2] (L) for k = 1, . . . , n can be determined by the first n—1 
functions F r ^ r /2] (L) for r = 2, . . . , n. 

Proof. The function Fl(0, 1, z) has the expansion, 

2n 

F L (0,l 7 z) = J2ML)z 2n - r , 

r=0 

where f r (L) are defined by f r = (— l) r / 2 F r r / 2 {L) if r is even, and f r = 2(— l)^ r / 2 ^F r r r / 2 i (L) if r is 
odd. From Lemma |2.I[ we have 

F L (0,l,z)=(z n + g 2 z n - 2 + ...g n ) 2 . 

where gu, k = 2, . . . , n are calculated from the Pfaffian of the matrix M{y, z). This implies that all 
g^s are determined by F r ^ r / 2 ] for r = 2, . . . , n. Then the next n functions F r ^ r / 2 ] for r = n+l, . . . , 2n 
are expressed by the quadratic forms of gk's, that is, F n+k ^ n+k y 2 ^ = ^k(F 2 ,i, . . . ,f n .[n/2]) with a 
polynomial function for k = 1, . . . , n. □ 
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Example 2.1. For the case of n = 2, we have the following two relations, 

^3,i = 0, F4 j2 = T-Ff,i ■ 

For n — 3, we have three relations: 

F4,2 — 7-^2,1) ^5,2=2^2,1-^3,1, ^6,3 = ~ -Pf,l • 

For n — 4, we have four relations, 

-^5,2 = 2^2, 1-^3,1: -^6,3 = — -Ffl + F2,iGa.2, -FV,3 = -^3,1^4,2, -^8,4 = G42 , 

where G 4 , 2 := |(fk,a - i^i)- 

Proposition 12 . 1 1 indicates that the curve C(L) has singularities of n double points at (Q,y, z), and 
that for generic L these are the only singularities, in which case the genus of C(L) is given by the 
Clebsch formula, 

(2.3) g(C) = i(2n - l)(2n - 2) - n = 2n 2 - An + 1 . 

Let us now find Casimirs for the coadjoint action Ad* for g E g: Recall that Ad* factors into an 
Sp-action and a G{-action. We first have: 

Proposition 2.2. The coadjoint orbit generated by Ad* P of P £ Sp(n, R) through X S sp(n)* is 
given by 

Ox = {Y esp(n)* : G fe (A) = C k (Y), k= l,...,n} , 

where Ck(X) is defined by the coefficients of Fx (2, 1, z) ; 

n 

F x (2, 1, z) = z 2n + Ck{X) z 2 <"- fe ) . 
fc=i 

Proof. With the Killing form on sp(n), one can identify sp(n)* with sp(n). An isomorphism of 
the identification may be given by 

X e sp(n)* = X' := X + JX T J e Sp(n) . 

Also with this identification, the coadjoint orbit Ox = {Aip(A) : P € Sp(n,R)} is isomorphic to 
the adjoint orbit O x , = {Ad P (X') = PX'p- 1 : P e Sp(n,M)}. Each point Y on the orbit O x can 
be then described by 

det(A + JX T J - zl) = dct(T + JY T J - zl) , 

which is equivalent to Fx{% 1, z) = F Y {2, 1, z). Note here that F x (2, 1, z) = F x {2, 1, -z). □ 
Thus {Ck{L) ■ k — 1, . . . , n} forms the set of Sp- invariant Casimirs, and they are given by 

k 

(2-4) C k (L)=Y / ^F 2k ^(L). 

j=i 

Remark 2.2. There is a corresponding proposition for the G{-orbit. In this case the invariants are 
generated by = det(L — zl)^) where (A)( k ) is the fc'th lower-left chop of A; the matrix found 
by deleting the last k rows and first k columns of A. The ratio of the first two coefficients of F^ for 
1 < k < n— 1 determines the Gj-orbit. However, for L lower Hessenberg, we obtain F^ = JXj=x a jt 
and Frty = , for k > 1, so that the Gt-invariants are trivial. 
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Thus we have n relations and n Casimirs in the set of {F ri k(L)} consisting of n 2 + 2n~ 1 functions. 
The total number of Hamiltonians for the Pfaff lattice is then given by 

n 2 -l. 

Recall that the L has 2n 2 + n — 2 free variables, and there are n Casimirs. This means that the 
Pfaff lattice (|1.6[) is a Hamiltonian system of degree n 2 — 1, and our goal is to give a complete set 
of n 2 — 1 Hamiltonians in involution, and to construct the angle conjugate to those Hamiltonians 
(action integrals). 

We now discuss the commutativity of the integrals, F r ^{L). Let us first note: 

Lemma 2.1. Let F, G be functions on q* . Suppose F is Gt-invariant, and G is Sp-invariant. Then 
{F,G} .(L) = O. 

Proof. We obtain 

{F, G} r (L) = (L, [n t VF, 7r e VG] - [tt sp VF, tt s9 VG\) 
= (L, [VF, tt e VG] - [^ 5P VF, VG]> 
= 0. 

Here the first term is zero by the G{-invariance of F, and the second is zero by the Sp-invariance of 
G. □ 

As an immediate consequence of this Lemma, we have 
Corollary 2.2. 

{F L (l,0,z),F L (l,y',z')} r (L) = 0. 

Proof. Note just that Fl(1,0,z) = det(L — zl) is both G(- and Sp-invariant, and FL(l,y' , z') is 
Sp-invariant. □ 

This implies {F rfi , F r > , k ,} g ,(L) = for r,r' = 2,3, ...,2n and k' = 0, 1, . . . , [r'/2]. In order 
to show the commutations {i*r,ki-^r',k'}fl* = for nonzero k,k' ', one needs to compute the case 
{Fl(1, y, z), Fl(1, y', z')} s * for nonzero y,y'. Let us consider the function for (l,y, z) £ C(L), 

H(y,z) :=hxF L (l,y,z). 

Using Vlndet(L) = L _1 , we have, by the chain rule, 

(2.5) VH(y, z) = (l- y)L{ Vl z)' 1 + yJL{ Vl z)~ T J , 

where we denote 

L(y, z) := (1 - y)L + yJL.- T J - zl = L y - zl 
and L~ T = (L _1 ) T . Now we can prove 

Lemma 2.2. For nonzero values of y — y' , we have 

{H(y,z),H(y',z')} 3 ,(L) = 0. 
Proof. Since H is Sp-invariant, we have 

{H(y, z),H(y', z')} r (L) = (L, [V% z), VH{y', z')]). 
We also have the following formula with L y := (1 — y)L + yJL T J, 

L = yL y > - y'L y 

y-y' 

From (|2~5l and [L y , L(y, z)] = 0, we have {H(y, z), H(y', z')} a * = 0. □ 
Corollary 12.21 and Lemma l2~2l now lead to 

Proposition 2.3. The integrals F r ,k all commute, i.e. {F r ,k, -Fr'.fc'jg* [L) = for all r, r' , k and kl . 
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One can then choose the following set of integrals as the Hamiltonians for the Pfaff lattice: 
(2.6) H(L) := { FrA L) : * Z ^ .' ^ } ■ 

The total number of Hamiltonians is n 2 — 1. In the next section, we construct the angles conjugate 
to those Hamiltonians. Then the independence of the Hamiltonians is shown by that of the angles. 

2.2. The angle variables conjugate to F r ^{L). Before we define the angle variables conjugate 
to the integrals F r ^, we start with a preliminary Lemma, and introduce several notions on the curve 
C(L): Recall the notation L(y, z) = (x — y)L + y.JL T J — zl, i.e. y, z) = det[L(y, z)}. 

Lemma 2.3. For generic C{L), we have dimc[ker(L(y, z))\ = 1, and the kernel is given by any 
nonzero row of the cofactor matrix L{y, z) c . 

Proof. For (1, y, z) G C(L), we have VFl(1, y, z) = Vdet[L(y, z)\ — L(y, z) c , from which we see 

F z (y,z) := £-F L (l,y,z) = -tr (L(y,z) c ) , 

and 

F y (y, z) := ^F L {l,y, z) = -tr (L{y, zf (JL T J - L)) . 

If L(y, z) c is the zero matrix then F y = F z — 0, which, for generic L, only happens at x = 0. If the 
cofactor of a matrix is non-zero, then the dimension of the kernel of the matrix is less than or equal 
to 1. However as the determinant of the L(y,z) is zero on C(L), the kernel has dimension strictly 
bigger than zero. Therefore the dimension of the kernel is 1. 

One may then show that if the dimension of the kernel of a matrix is one, then any non zero 
row of the cofactor matrix is a basis for the kernel. This is a straightforward computation using 
Cramer's rule. □ 

We define the divisor T>(L) on C(L) such that with the eigenvector map / : C(L) — ► CP 2n_1 ;p 
f(p), with L(y,z)f(p) = (L y — zl)f{p) = 0, the first element of f(p) vanishes, denoted as fi(p) = 
(ei,f(p)) = for p e 'D(L). In terms of the cofactor L(y,z) c , the V{L) can be expressed as 
V{L) U lD{L) = {(L(y, z) c )i ; i = 0}, where iD(L) is the image of V{L) under the involution t : 
(x, y, z) — > (x, x — y, —z). For some yo ^ 0, 1 wc denote Pi — (1, y^, Zi) G C{L), i = 1, 2, . . . , In (the 
fiber of C(L) over y ) and Qi = bP t = (1, 1 - y , -z{) for i = 1, . . . , In, i.e. L yo f(Pi) = z i f(P l ) and 
L\- yo f(Qi) = —Zif(Qi). Then we consider the dense set of matrices L satisfying 

(a) V(L)n{F z (y,z) = 0} = <&, 

(b) T>(L) is a finite set of C(L), 

(c) V(L)n{{(L(y,zy) 1 . 1 } z =0} = ®. 

(d) V{L) n {y = 0} = 

(e) {P l7 Qi} n V{L) = for (1, y , z ) $ C(L), 

(f) {Pi, Q t } n {F z {y, z) = 0} - for (1, y , z Q ) £ C{L). 

(g) The only singular points of C(L) are at x = 0. 

Here [(L(y,z) c )i : i] z = d[(L(y, z) c )i t i]/dz. The condition (c) gives the genericity of the divisor 
T>(L). By the Riemann-Hurwitz theorem, we also note that the degree of the divisor T>(L), d(D) := 
degV(L), is given by 

d(V) = g(C) + 2n - 1 = 2n(n - 1). 
Now we consider a family of differentials on C(L) which define the angle variables conjugate to 
the action integral F r ^{L), 

rt7 > , T y. < x dy Jfc = 0,l,...,n-2, 

(2J) ^^=^^F^ry fOT V- = 2fc + 2,...,2n. 
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where g r ,k will be specified later. By the Riemann-Roch theorem, there exist g(C) = 2n 2 — An + 1 
holomorphic differentials on the curve C(L). As will be shown below, the set of those differentials 
consists of n(n — 2) odd differentials and (n — l) 2 even differentials with respect to the involution 
l : (1, y, z) — > (1,1 — y, —z). We will use the odd holomorphic differentials to give the angle variables 
conjugate to F ri k for k = 1, . . . , n — 2 and r = 2k + 2, . . . , 2n. For these holomorphic differentials, 
we use polynomial functions of y and z for g r .k, and since deg(F z ) = 2n — 1, 

deg(g r ,k) < 2n - 3. 

The angles generated by the odd holomorphic differentials account for only n(n — 2) of the n 2 — 1 
necessary angles. The remaining 2n — 1 angles will be generated by meromorphic differentials uo r fi 
for r = 2, . . . , 2n with poles at y — 0. This meromorphicity at y = 0, together with the involution 
on C, explains why we need condition (d) and yo ^ 0,1. 
Let us denote the singular locus of C(L) by 

P oo (L) = {(0,l,7?i) : 1 < i <n} , 

and denote 

Vi{L) = (l.Oi.fti) S 2?(£). 

Note that T> OQ {L) is invariant under the coadjoint action of Sp(n), as the singular points are a feature 
of the curve C(L), while T>i{L) evolves. Then we consider the integral of u> r k from VaoiV) to T>(L): 

d(p) Vz 

i=l 

This defines a map, C(L) d ^ — * Jac(T>) where C d( ^ l; '' ) is the d(P)-symmetric product of C(£) and 
Jac(T>) associated with the divisor T>(L) is an extended Jacobian defined below (see also [8]). Then 
we show that <& r ^(L) gives the canonical angle conjugate to the integral F rj k(L), that is, we determine 
g r ,k(y,z) so that we have 

Proposition 2.4. For any F r > # E H(L), we have 

(2.8) {®r,k,Fr>,k>} s , (L) = 6r,r'Sk,k> for 

The total number of those Q r k is n 2 — 2n, and the other 2n — 1 angles given in Proposition 12.51 
below are conjugate to the Chevalley invariants F r- o for r = 2, 3, . . . , 2n. 
To show Proposition ^. 41 we first note: 

Lemma 2.4. Let H(yo, zq) = ]nFz,(l,yo, Zo) for (l,yo, zq) ^ C(L) and yo ^ 0, 1. Then we have 

(2.9) {$ r . fc , H(y , z Q )} (L) = — — (G r .k(y , z ) + G r . k {l - y a , -z )) , 

F L {l,y ,z a ) 

where G r ,k(yo, zo) = yo(l - yo)(2y - l)g r ,k(yo, zo). 

Proof. The commutator can be written as the flow generated by the Hamiltonian H(yo, zq), 

(2.10) {$ r , fc , H (y , zq)} s . (L) = j t ®r,k{L) 

(2n) - y 9r,k(ai,0i) daj 

' ' f-f Fz(ai,bi) dt ' 



$ r ,fe(L) = / W r ,fc(L) = 



fc = l,...,n-2, 
r = 2& + 2,...,2ra. 
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Recall that on T>i — (l,Oj,&» = z(ai)), the eigenvector f(T>i) for L ai f(T>i) — bif(T>i) satisfies 
(ei,f(T>i)) = 0, i.e. T>i <E T>(L). The evolution of f(T>i(t),t) is given by 

(2.12) ^f(p h t) = -B y f(p i ,t) and f(V i ) = f(V h 0). 

Here the generator B y can be found from the Lax equation = [L y , B y ], which is the isospectral 
condition for L y f = zf: Recall that the Pfaff lattice with the Hamiltonian H := H(y ,z ) is 
^ = [L,7r f Vff ] (see ([O])). Then we have 

(2.13) ^ = (l-y)[L,7r e VHo]+yJ[L 1 Tr t VH ] T J 

at 

(2.14) = [L y , ^ { VH ] + y[JL T J, JineVHofj - ^ { VH ] • 

We now compute JL T J and J(tt^W Ho) T J — tt^Hq: Let us first note by the chain rule, 

ViJ = (1 - yo)L(yo, zoy 1 +yoJL(y ,z y T J . 

Then we have 

TreVffo = (1 - 2y )7r e (i(j/o, Zo) _1 ) , 

which leads to 

(2.15) J(tt ( VH ) t J - 7r e Vffo = (2yo - 1) (L(y , zq)" 1 - JL(j/ , z y T J) . 

We also note from L y = (1 — t/)L + y.JL T J and the same with y — yo that JL T J can be expressed 
as 

(2.16) JL T J = L a — Ly,, = — -L y — -JLE J . 

y-yo 1-2/0 2/ -(1-2/0) tf-(l-ito) yo 

Substituting ()2.15p and (|2.16| into (|2. 13|) . we have the Lax equation 

dL y 
~dt 

with 

B y = ir ( VH + (2y - 1)2/ (^—^L(y , zq)" 1 ^ r JL(y , z Q y T J 

Now from ([2~T2"j) and (ei, f{V l (t), t)) = for each t, we have (e x , §£ + §£:^) = which leads to 



~d/~ 



{e u B v f<Pt)) 



t=o [ei,d y f(Pi)) 

Since 7rjViJo is a lower triangular matrix, (ei,TTtVHof(T>i)) = 0. Then the RHS of (|2.10[) becomes 
(2.17) 

9r,k( a i' b i) ( 2 2/o-l) / faj(l-2/o) r / N-i owo 



T7 TT7 5 j/-n ^ ei ' 1 L{yo,z ) -JL{y ,z ) J > f[T>i) 

. F z (a h bi) {ei,d a J{T>i)) \ [ a, ~ y m - (1 - y ) J 

Recall that (ei, f(T>i)) = for T>i = (1, a i: b{) £ T>(L). Also we have that L(y , z )f(yo,Zj) = (z 



zo)f(yo,Zj) and JL(y , z ) T J f(l-y , -Zj) = L(l-y ,-z )f (l-y ,-z j ) = -(zj- z )f(l-y , -Zj) 
for j = 1, . . . , 2n. Then using the residue theorem, equation (|2.17[) can be written by 



-1/0(1 -K>) (22/0-1)^ 



2n 

9r,k(yo,Zj) 1 . ffr,*(l -2/0,-Zj) 1 



^(2/o, % - z F z (l - yo, 2j - z 
Again by the residue theorem, we finally obtain 

Tit \\ 2/o(l ~ 2/o) (22/o - 1) , 

{$ r , fe ,ff(2/ ,^0)} n » = FTTi S (ffr,fc(2/0,20)-.gr,fc(l"2/0,-^0)) , 
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which is ([279]) . □ 

We now prove Proposition ^. 41 

Proof. With (|2.7p . we choose the following polynomials for g Ti k(yo, zq): 
For the case r = 2s + 1, we choose, 

I \ 1 / wfe-l 2(n-s)-l r I K = 1, . . . , n — 2, 

32 S +i,fc(2/o,2 ) = o(yo(l - yo)j V for < 

^ I S = K + 1, . . . ,71 — 1. 

Then we have 

(2.18) {*a^i,*,tf(lfo, «>)}*•(£) - ~p~FT~ y 2s+hk) (l,yo,zo). 

F L {l,yo,z ) 

For the case r — 2s, we choose 

/ \ (yo(l - 2/o)) fe_1 ,„ \\a-k\ Hn-s) f Jfe = l,...,n-2, 

92sAyo,z ) = — ^ 7\ — (i - (%>(i - yo)) )v for i , , -, 

2{2y -l) [s = k+l,...,n, 

which gives 

(2.19) {$ 2Sik ,H(y ,z )} B *(L) = — -i (> 2s < fc ) (1, y , z„) - 4 r -y 3 '-«)(l, »>,*>)) • 

fill, 2/o, zo) v ' 

On the other hand, using (|2.2p for Fl(1,j/o, z o), we have 

-. 2n [«-'/2] 

(2.20) {$ r , fc , J ff(2/o,zo)} fl .(i) = ttt^ r E E /'"''(^o^ol^^r'^VW- 

Comparing this to (|2.18[) and (|2.19[) , we get the results of Proposition ^. 41 Notice that F2 S , S is not a 
Hamiltonian in our list, hence the presence of the <^( 2s ' s ) term in (|2. 19[) does not affect the statement 
of the Proposition. □ 



Remark 2.3. The function <?2s, fe is indeed a polynomial and when inserted into (12. 7|) gives u>2r,k a 
holomorphic one form for 1 < k < n — 2 and k + 1 < s < n. 

In order to find the angles $ ri o(F) conjugate to F r ^(L), we consider the following integral of the 
meromorphic differentials u rj o(L) with a pole at y = 0, 

$,o(L) = f W ^ o( L), ^o(L) = MMl^vfalf) . 

j (L) y F z {y,z) 

where m r (y, z) and n r (y, z) are polynomials and m r (0, z) =^ 0. Then one can show: 

Proposition 2.5. For F r ' : fc' € Ti-(L), we have 

{$r,o,-Fr',k'} B * ( L ) = tir,r'8k>,o IOT r = 2, 3, . . . , 2n. 

Proof. As in the previous case, we compute {$r,o, H(yo, zo)} g * (L) with H(yo, zq) — ItiFl(1, yo, zq) 
for (l,yo,zo) 4- £-(£)■ Then following the previous calculation, we obtain 

{*r.o, H(y , z )} (L) = — ^ — l — [{(1 - y )m r (y Q , z ) - y m r (l - yo, -z )} 
8 F L {l,y ,z Q ) 

+2/o(l - Vo) {n r (yo, zq) - n r (l - y , -Zq)}] . 
Now we choose the following m r (yo, Zo) and n r (yo, zq): For r = 2s + 1, we choose 

( I \ 2(n-s)-l 

m 2s+1 ( yo ,z )=z , for s=1) n _ L 

n 2s+ i(y , ^o) = 0, 
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Then we have 

(2.21) {*2 5 +i,o, H(y , z )} r (L) = —± ^ 2s+1 ' 0) (1, y , z ) 

For r — 2s, we find that the polynomials 

m 2s (yo, Zq) = (2y a - l)z^ n ~ 

( \ ^ , ,. n \\s~l\ 2(n-s) for s = 1, . . . ,n, 

n 2s {yo, z ) = (1 - (4y (l - Vo)) ) V 



give 



(2.22) {<Z> 2Si0 ,H(y ,z )} g *(L) = — -i U*°M (l, yo , z ) - 4 y 2 ^(l, y ,z )) . 

Now comparing (|2.20|1 for {$ r ,o, H{y , zq)} s » (L), we obtain the desired results as in the previous 
case. □ 

Propositions 12.41 and 12.51 provide a proof of the main theorem in this section: 



Theorem 2.4. The Pfaff lattice fl.3\) with L 6 Zr is a completely integrable Hamiltonian system 
with n 2 — 1 Hamiltonians given in \2. 6]) and n Casimir functions of 

2.3. Extended Jacobian. Here we just mention that the Pfaff flows are linearized on the extended 
Jacobian defined below. First recall that the curve C(L) has the genus g(C) = 2n 2 — 4n + 1 (recall 
(|2.3|) ), and the positive divisor T>(L) has degree d(T>) = 2n(n — 1). There are g(C) number of 
holomorphic differentials (the Riemann-Roch theorem). We have n 2 — 2n differentials ujr^ defined in 
(|2.7j) for k = 1, . . . , n — 2 and r = 2k + 2, . . . , 2n, which are odd forms with respect to the involution 
i : (1, y, z) — > (1,1 — y, —z). In addition to those, we have other (n — l) 2 holomorphic differentials 
uifi k , which are even forms defined as follows: For r' = 2s + 1, 

which are just w^s+i fc = i^V — l) w 2s+i,fc) and the total number of these is ("~ 1 K"~ 2 ) _ p or r ' _ 2s, 

< fe = (y(l-y)) & -^ 2( "- s) I 



: l,...,n- 1, 
fc + 1, . . . ,n, 



which give "*-" 2 differentials. Then we have an extended Jacobian Jac(T>) associated with the 
divisor T>(L) whose dimension is d(T>), 

(2.23) Ja~c(V) = t $ r<k (L) = / w r , fc (L), * r ',fc'(£) = / wi 1 *)) ■ 

Note here that the differentials u> r> k(L) also include the meromorphic one forms defined for k = 
and r — 2,3, ... , 2n. 

Remark 2.5. Since the & r ',k' m (|2.23p are even with respect to the involution t, they are constant 
under the flows of every Hamiltonian F rk . 



Then the commutation relations in Propositions 12.41 and 12.51 mean that the flows of the Pfaff 
hierarchy are linear on the Jacobian. 

The extended Jacobian Jac(T>) of genus g = g(C) is well defined up to the choice of paths in 
the integrals in (|2.23[) : the path of integration may be deformed by an element of the homology 
Hi(C \ {y — 0} , Z). Choose a basis of H\(C \ {y = 0} , Z) in the following way: Let {A,-, B{\ 1<i<g 
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Figure 1. A real section of C(L) with n = 3 together with the real points of T>(L) 
for f > 0. In this case the six eigenvalues of L are {—1,-0.5,-0.25,0,0.5,1.25}. 
The right figure is the part of the left one near the fiber y = 0. In particular, note 
that the rate at which the points in T>(L) approach the eigenvalues depends on the 
size of the corresponding eigenvalue. 



be a canonical basis for i/i(C,Z) with intersection product Ai ■ Bj — Sij and let {Cj} 1 < J < 2 „ be 
(homologically trivial) cycles about the points (1,0, Zj) with a positive orientation, where the Zj are 
the eigenvalues of L. Then the Jacobian map (|2.23p is well defined up to the lattice T generated by 
the columns of the (g + 2n — 1) x (2g + 2n) period matrix 



n 



Pi o 

P 2 N 



where the rows of II correspond to the differential forms, and the columns of II correspond to the 
cycles. Pi is the g x 2g period matrix of the holomorphic one forms with respect to the canonical 
basis {A i: B{\. P 2 is the (2n — 1) x 2g matrix of periods of the meromorphic one forms uj r ^ over the 
{Ai, Bi}. While TV is the (2n — 1) x 2n matrix where the entry 



N 2s -x,j - 



2(n-s) 



3,(1, 0,Zi) 

is the residue of u>2s,o at (1, 0, Zj) for s = 1, . . . , n, and the entry 

2(n-s)-l 
Z ■ 

N 2s i = 2m^- - 

is the residue of W2 S +i,o at (1, 0, zj) for s = 1, . . . , n — 1. Note that N = D1VD2, where D\ and D 2 
are diagonal matrices, and V is a (2n — 1) x 2n Vandermonde matrix with the entries (V)ij = z 1 ^ 1 
for 1 < i < 2n — 1 and 1 < j < 2n. The residue theorem and a standard Vandermonde argument 
imply that ker(iV) = span c { (1, 1, . . . , 1) }. Therefore the rank of r is 2g + 2n — 1. As a real space, 
C 9+2n_1 /r has 2g + 2n — 1 compact dimensions and 2n — 1 noncompact dimensions. 

If the eigenvalues Zj are all real, then the entries of iV are pure imaginary. One then checks that 
d j['° = {&r,o, Fr,o} = 1 so that $ rj o evolves only in the real direction. Therefore as N has imaginary 
entries and full rank, and Pi generates a rank 2g lattice, the evolution of $ r .o is entirely in a copy 
of M rather than S 1 . This implies that the Chevalley flow induces a noncompact torus action. 

Numerical calculations with n = 2 and 3 show that the divisor T>(L) contains some of the points 
(1,0, Zj), as ti — > ±00 (see Figure Q}. One then notes that the meromorphic one-forms used to 
produce the angles <& ri o conjugate to F r $ posses poles at (1, 0, Zj). The 2n— 1 noncompact dimensions 
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of the Jacobian produce the fixed points of the noncompact torus actions generated by the Pfaff 
lattice hierarchy. 

One can construct the solutions of the Pfaff lattice using the Riemann theta function through the 
Jacobi inversion problem. However we will not do this here. Instead we use the matrix factorization 
to construct the solutions in the next section, and will describe the (real) isospectral variety using 
the moment map for the noncompact torus actions in Section [5l 



3. Matrix factorization and the t-functions 



In this section we derive the solutions of the Pfaff lattice hierarchy in terms of the SR- factorization 
of a matrix in SL{2n). We define the r-functions whose derivatives generate solutions. In Section 
E] we use the r-functions to compute the fixed points of the Pfaff lattice by computing asymptotics 
as t — > ±00. We conclude this section by examining the r-functions in the GOE and GSE-Pfaff 
lattices. 

As in the case of the Toda lattice, we consider the following matrix factorization for e e ^' L ^ € 
SL(2n,R) with the initial Lax matrix L(0) and 9(t,L(0)) := Y^lJ * tjL(0) j , 



(3.1) 

where P £ Sp(n, I 
G t = 



J(t M 0)) = Q( t )-lp( t ) , 



and Q G Gf. The groups Gj and Sp(n,M.) are defined by 



fhxh 2 



2 \ 
2 



5p(n,R) = {P e SL(2r 



h n h) 
:PJP T 



fc=i 



dimGe = 2n 



J} , dim Sp(n, R) = 2n 2 + n . 



This factorization is called the SR factorization (a symplectic version of the QR factorization), and 
the set of matrices having an SR factorization is dense in SL(2n, R), but not in SL(2n, C) (Theorems 
3.7 and 3.8 in [1]). A standard computation starting from (|3.ip shows that the matrices P(t) and 
Q{t) satisfy 



(3.2) 



dP (Ti\-P 9 ® 

^7 = 7rsp(i)P ' dT 3 



-ir t (L J )Q, j = l,. 



,2n-l. 



and Q(0) = P(0) = I, the identity matrix. The evolution of the Pfaff matrix L(t) is then given by 

(3.3) L(t) = Ad Q{t) L(0) = Q{t)L{Q)Q{t)- x . 

One should note that the factorization (13. ip is not always possible, and in general, for some t, we 
have 



wQ 1 P for some w G S2r. 



with the symmetric group S% n . For example, consider the case L(0) = Co, the companion matrix 
(|1.8p with 7 = 0, and set t% = 1, = |, ts = |, which then gives 



(3.4) 



J(t,L(0)) _ 



(1 


1 


1 


2 ^ 





1 


1 


1 








1 


1 










V 
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A factorization of this matrix is given by 



tf'C 



1 P 



A) 





1 














v\ 




(° 





1 


1\ 











1 







1 



















1 


i 













1 


1 


1 







1 


1 








V> 


1 









vl 














1 





0/ 



Here w is the permutation (1, 3)(2, 4). Note that there are other choices of w, for which w~ 1 e e ( t ' L (°^ 
has other SR- factorizations. The factorization e e ^' L (°^ = wQ~ 1 P with w ^ id corresponds to a 
singular solution of the Pfaff lattice, and will be discussed in a separate publication (see [6] for the 
case of the Toda lattice where the factorization is given by the Bruhat decomposition). 

3.1. Moment matrix and the r-functions. Let us define a 2n x 2n skew-symmetric matrix 
(3.5) M(t) = e e ^ L ^Je e ^ L ^ T =Q- 1 JQ- T for Q £ G t , 

which we call the moment matrix of the Pfaff lattice of (|1.6p . This Cholesky-type (or skew Borel) 
decomposition (|3.5[) represents the Gauss elimination of a skew symmetric matrix M. Again one 
should note that in general the decomposition has the form, for some t, 

M(t) = wQ~ 1 JQ~ T w T for some w £ S 2n ■ 

In the above example (|3.4p , the moment matrix has the form, 

/0 -1 l\ 



M = 







V 



i 
i 



(In this paper, we leave blank the lower triangular part of skew-symmetric matrix.) Notice that 
the element m^a = 0, and this is the reason why one needs a permutation w ^ id to have the 
decomposition, M — wQ~ 1 JQ~ T w T with Q £ Gt. 

Let us normalize Q by a block diagonal matrix H :— diag 2 (/ii/2j ■ ■ ■ , h n l2) with Ofc=i = 1 and 
Q = R- X Q so that 

(h 2 ■•• 2 ^ 
* h ■ ■ ■ 2 



= hn + R , 



where R is a 2 x 2 block lower triangular matrix. This corresponds to the group factorization, 
G( = A ■ Gt, where A is the set of diagonal matrices of the form of H and Gt is the set of 
matrices of the form of Q. The Gt is a unipotent subgroup of SL(2n, K), and can be identified 
as the top cell (open dense subset) of the flag variety SL(2n,M.) / 'P 2 x2 where P 2x2 is the parabolic 
subgroup of 2 x 2 block upper triangular matrices. This implies that the canonical morphism 
if> : SL(2n,M) — > 5L(2n,M)/P 2x2 gives an open immersion for the restriction of ip on Gt, i.e. 



(3.6) V : Gt £«i , 

where Ej^ is the top cell of the flag SL(2 Tl 5 M ) / -T 2 x 2 • 
With the definition of Q, i.e. Q = HQ, we have 



Lemma 3.1. The matrix Q satisfies 



dQ 



— — R ^ 
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where Bj is a 2 x 2 block lower triangular matrix given by 

B J =H{(V)--J(LJ)IJ) H-K 

Also the matrix H satisfies 

^=\{{IJ) -J{U) T J)H, 

Proof. From (|3.2p . we have 

Since Q^ 1 ^- € 0- and ^-H' 1 £ q q , calculating the projection nt(lJ) the statement follows. □ 

The elements hk in H can be calculated from the Pfaffians of the matrix M: Let M 2 fc be the 
2k x 2k principal part of M, the 2k x 2k upper-left block of M. We then define the r-function Tik 
as the Pfaffian of M2k, 

fc 

(3.7) r 2k = pf(M afc ) = P f((i/ 2 J) 2fc ) = J]^, fc = l,...,n. 

i=i 

where (-ff 2 J) 2 fc = diag 2 (/i 2 J 2 , . . . , ft. 2 J 2 ). The explicit form of the Pfaffian of a matrix in (|3.7[) is 
given in (|3.15|) . 

We wish to describe the fixed points of the Pfaff lattice by computing the limits of the flow as 
tj — > ±oo; however for the Pfaff variable L € ^(7) the flows are unbounded in those limits (the 
Chevalley flows correspond to non-compact torus actions on the extended Jacobian, see Section HO)) . 
We then normalize these flows by conjugating L by H: define L as 

L = HLH- 1 = QL(0)Q-\ 

from which the elements ak in the super diagonal of L are given by 

/q o\ ((\\ hk + l ln\ V /T 2'c+2 r 2fc-2 

(3.8) a k = a k (0)— — = a k (0)- , k = 1, . . . , n - 1 . 

«fc T 2fe 

All other entries of L can also be described by the set of r 2 fc functions as will be shown in Section |4] 
(see also [3]). In this paper, we consider the case where Ofe(0) are nonzero, and we normalize L by a 
conjugation so that the elements on the super diagonal are all 1. A Pfaff lattice having a/-(0) = for 
some k splits into smaller lattices (sub-systems) which appear as the boundaries of the isospectral 
variety of the generic flow. We then define the isospectral variety of the normalized Pfaff lattice 
with matrix variable L as 

(3.9) Z R (7) = {lg4(t) : ^fc,fc+i = 1, k = 1, . . . , n - l} . 

Now the asymptotic matrices of L(t) as tj — > ±00 are bounded, lie in Z^j), and are fixed points of 
the Pfaff lattice. 

Remark 3.1. In the case of tridiagonal symmetric Toda lattice, the symmetric moment matrix and 
r-functions are defined as follows: Let Lr{t) be an n x n Lax matrix for the Toda lattice of sl(n, K). 
We have the QR factorization, 

e 9(t,L T (0)) = r (t)~\(t) r(t) E B+, q(t) E 0(n) , 

where B + is the Borel subgroup of upper triangular matrices. The (symmetric) moment matrix My 
is then defined by the Cholesky factorization, 

M T = e 2e ^ LT ^ = (r- 1 q)( (Z T r- T ) - r" 1 ^ . 
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The r-functions are defined by 

(3.10) rl oda = det((M T ) fe ) , 

where {Mt)r is the k x k upper left submatrix of the Mr matrix. We also obtain similar formulae 
to (|3.8p for the off diagonal elements in Lt matrix, which are also denoted as at (see e.g. [7]). 

We note that T 2 k can be parametrized by a constant skew-symmetric matrix determined by the 
initial matrix L(0), which we call the B-matrix, as follows: With the normalization etfc(0) = 1 (see 
above), the elements in the super-diagonal of L(0) are all 1. One can then write an eigenmatrix $o 
for L(0), in the form $o = ^aV where € G{ with diag^o) = and the Vandermonde matrix V 
is given by 

1 



V 



( 1 



1 \ 



Z2 



,2n-l 



Z-2r. 



_2n-l 
z 2n 



Then we have 

M = Soe^'^o^o^'^o 
= y E(t,A)BE(t,Afy% . 

Here the £>-matrix is defined by 

(3.11) i? = $- 1 J$o T , 

A =diag(zi, . . . , z 2n ) (i.e. L(0)$o = 3>oA with the eigenvalues Zi), and E(t, A) 



y e 0(t,A) = e 6(t,c-,)y 



is given by the Wronskian matrix of {Ei(t), . 

( E x 

(3.12) E(t,A) = 



. , E2 n (t)} with respect to the t\ variable, i.e. 



E 2 
EL 



E-2n 
E 2n 



(2n-l) 



E. 



(2n-l) 



E. 



(2n-l) 
2n 



with Ek(t) = expd(t, Zk) for k = 1,. 
£'(<:, A) satisfies the linear equations, 



, 2n, and 



0') _ &E k _ 



z° k Ek- Note also that the matrix 



(3.13) 



J-£(i,A) = C*£(i,A), 



E(0,A) = V, 



where we have used VA = C 7 V with the companion matrix C 7 . Since 'J'o does not affect the value 
of r-functions, it is convenient to define M := \E , ^" 1 M\I/ ( ^ T = EBE T . The M then satisfies 



(3.14) 



^ = C*M + M{C*) T , M(0) = VBV 1 
Otk ' 1 



This equation can be considered as a linearization of the Pfaff lattice on the space of skew symmetric 
invertible matrices as described in Section O The B- matrix then gives the initial point of the flow, 
hence it determines the symplectic leaves of the foliation with the Hamiltonians, F r ^k{L). 

We remark here that the Pfaff flow can be linearized on SL(2n,M.)/ Sp(n,M.) as an image of the 
companion embedding, c 7 : Z^(j) — ► SL(2n, M.)/Sp(n, R), that is, we have the commuting diagram, 

£(0) Cy » mod Sp(n,R) 



L(t) e 0(*,c 7 )^-i m0 dSp{n,R) 
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Here we have used L(0) = ^oC^q 1 and L(t) = Ad Q{t) L{Q) with e e ^ L( ^ = Q(t)- 1 P(<). Note here 
that the Pfaff flow on SL(2n, M.)/Sp(n, R) satisfies the same equation as E(t,A) in (|3.13l) . 
Now we can express the r-functions in the form with 2k x 2k upper-left block of M, 

(3.15) r 2k = pf(M 2 fe) = ^2 • • ■, i k,jk)m il ,j 1 m i 2,j2 ' ■■ m i k ,j k ■ 

where the sum is taken over the set I 2 k = {1 = «i < • ■ • < ife < 2k, i s < j s . s = 1, . . . . k}. The 
coefficient a(ii,ji, ■ ■ ■ > ik,jk) is a parity of the permutation, i.e. 

. (\ 2 ■■■ 2k -1 2k 
a := sign . . 

\*i Ji ■ ■ ■ «fc Jk 

For example, we have 

T 2 = TOl j2 , 

r 4 = toi j2 to 3i 4 - mi :3 m 2i 4 + wi,4"i 2: 3. 
The r-functions are also obtained as the exterior products of the bi- vector £! 2 := rriijei A ej, 

l<i<j<2n 

i.e. 

r 2fc = (A k n 2 ,e 1 A ■ ■ • A e 2fc ) , 

where (•, •) is the usual inner product on A 2fc R 2 ™. This indicates that the nonzero r-functions give 
the generic orbits of the Pfaff lattice, that is, those orbits are on the top cell of the homogeneous 
space SL(2n)/Sp(n) (see Section [5]). 

The entry m<j of the moment matrix M = EBE T is explicitly expressed as 



E 

Kk<K2n 



/ J ~k,l 



(3.16) 

" hkAz k ziT\z 3 r - 



\<k<K<2n 



From this expression, we have 

(3.17) 1,3 — m i+ k. 3 + mij + k k = l,2, ... 

atk 

This, of course, is the same as (|3.14[) . but notice that (|3.14[) includes the characteristic polynomial 

<hn{*) = Z 2n + Et n =2(-W n - k = 0. 

Then it can be shown that the r-functions satisfy the Hirota bilinear form [13j . 

(3.18) (-^D 1 D 3 + Di + 3Dl)T 2k -T 2k = 24r 2k _ 2 r 2k+2 , k = 1, 2, . . . , n - 1, 

with t = 1. Here D k is the Hirota derivative with respect to t k , D k f ■ g :— (d tk — d t i k ) f (t)g{t')\ t=t i . 
Note here that the first equation in (|3. 18|) with r 2k+2 T 2k ^i = is the usual KP equation, and the 
T 2k for the KP equation is given by the Wronskian form (for example, see [18:). One should note 
in (|3.18p that r 2k are all generated by r 2 . This is similar to the case of Toda lattice where the 
r-functions T^ oda satisfy 

n 2 Toda Toda _ Toda Toda . _ i i 

u \ T k T k — T k~l T k+1 ' K — 1, 1, 

with TQ° da = 1. From this recursive equation, one can show that rJ oda are given by the Wronskian 
(Hankel) determinants, 

T Toda = Wr(j Toda^ ^Today ^ ^ ^Toda) (fc-1) j ^ 
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where (r 1 " oda )^) = T^ oda f dt\. This expression of the r-function then agrees with (13 . 1 0[) in Remark 

ED 

Remark 3.2. The system f|3 . 1 8[) has been proposed as a first member of the DKP hierarchy in [T5] . 
In (T3] , Hirota and Ohta introduced the system as an extension of the KP equation and called it the 
coupled KP equation. The system (|3.18p was rediscovered as the Pfaff lattice describing the partition 
function of a skew-symmetric matrix models in [H [21 H7|. In particular, Adler et al in 2 discussed 
the lattice structure of the system and formulated it as the Pfaff lattice with the connection to the 
Toda lattice. The system (|3.18j) was also found as a charged BKP hierarchy describing an orbit of 
some infinite-dimensional Clifford group action in [16] . 

3.2. Foliation of the phase space by F r k(L). The integrals F r {L) (the Chevalley invariants) 
define the isospectral manifolds ^(7). We have shown that the ^(7) embeds into the space of 
skew-symmetric invertible matrices. This space is then foliated by the additional integrals F r k(L) 
for k > (found in Section^. Our point of view is that the i?-matrix of (|3.11[) gives coordinates 
on the space of skew-symmetric matrices. We can then write the constants F r ^(L) in terms of the 
eigenvalues {zi : i — 1, . . . , 2n} and the B- matrix: We have 

Proposition 3.1. 

F L (x, y, z) = — — det[(x - y)AB - yBA - zB\. 

Proof. Using L = QL^Q- 1 and M = Q- 1 JQ~ T , we have 
Fl(x, y, z) — det[(x — y)L + yJL T J — zi] 

= det[{x - y)QL{0)Q- 1 + y.JQ~ T L(0) T Q T J - zi] 
= det[(x - y)L(0) - yML(0) T M _1 - zi] . 
By factoring M out and setting t = 0, this can also be expressed as 

det(M(0)) F L (x, y, z) = det[(x - y)L(0)M(0) - yM{0)L(0) T - zM{0)} 
= dct[(x - y)AB - yBA - zB] det($ ) 2 • 

In the last step, we have used L(0)$ o = $oA, M(0) = $ q B$q T . Then noting det(M(0)) = 
det(-B) det(<I>o) 2 , we obtain the result. □ 

We denote pf(zi,?2, as the Pfafhan of the k x k skew symmetric matrix found by taking 

the i\, i 2l . . . , ik rows and columns of the £?-matrix, e.g. pf (B) — pf (1, 2, . . . , 2n). We note that 
pf(£) = [det(U)]" 1 = [U l<3 ( z 3 - Zi)}' 1 sm cc 1 = det(M(Q)) = det($ ) 2 det(B). 

Example 3.3. The case of n — 2: The F r ,h(L) are expressed as 

^2,0 = ZlZ 2 + ZlZ 3 + Z1Z4 + z 2 z 3 + z 2 z A + z 3 z±, 

F 2 ,i = 2 (61,263,4(21 + z 2 f - 61,362,4(21 + z 3 ) 2 + 61,462,3(21 + 2 4 ) 2 ) /pf(B), 
^3,0 = (21 + 2 2 ) (21 + 2 3 ) (zi + z 4 ), 
^3,1 = 0, 

F 4 , = Z1Z22324, 

^4,1 = - (61,263,4(21 + 2 2 ) 2 (Z1Z 2 + 23Z4) - 61,362,4(21 + Z 3 ) 2 (ziZ 3 + Z 2 Zi) 
+61,462,3(21 + 2 4 ) 2 (ziZ 4 -I- 2 2 2 3 )) /pf(-B), 

with F 4 .2 = \{F 2 .\) 2 (see Example 12. ip . Here bij are the elements of the B-matrix. Recall that 
F r _Q for r = 2,3,4 are the Chevalley invariants which do not depend on the i?-matrix, and C\ — 
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4-F2,o + ^2,1, C*2 = 16^4,0 + + F4.2 are the Casimirs. This gives F24 = £1 , F4.2 = £2 , and in 
addition we have pf(-B) = 61,263,4 — 61,362,4 + 61,462,3 — £3, where ci, £2, and £3 are constants. This 
system of equations is linear in 61,263,4, 61,362,4, and 61,462,3. For generic z\, Z2, and Z4 it is a 
non-singular system. The foliation of the isospectral variety Zg_(~/) is therefore parametrized by 

ci = 61,363,4, c 2 = 61,362,4, c 3 = 61,462,3 , 

where Ci, C2, and C3 are constants. Note that as the pf(£>) 7^ not all c^'s are zero. This foliation is 
related to the cell decomposition of the space of non-singular skew-symmetric 2n x 2n matrices. 

The case of n = 3: We find for instance that 

(3.19) F 2 ,i= (-l) i+J ' +1 &ijPf(l,2,... J,...,j,...,2n)(z l + z,) 2 / P f(i?), 

l<i<j<2n 

where i means that these elements are skipped in the list. Our hypothesis is that this formula is 
true for general n. 

Remark 3.4. One should compare these results with the formula computed in 9] ((23) in p. 197) for 
the fc-chop integrals of the full Kostant Toda lattice which are found to be rational functions of the 
Plucker coordinates ttj and 7r} with J C {1, . . . ,n} and \ J\ = k on the flag manifold SL(n, C)/B + 
with the Borel subgroup B + . 

3.3. Examples from the matrix models. The Pfaff lattice was introduced as an integrable lattice 
whose r-functions are the partition functions for the matrix models of the Gaussian orthogonal 
ensemble (GOE) and Gaussian symplectic ensemble (GSE) [T]. The partition function of a random 
matrix model is given by the matrix integral, 

(3.20) r(t u h,...)= [ e M-v(x)+«(t,x)) dX ; 

J Mn 

with Haar measure dX on the matrix ensemble M. jf, and some potential function V(X), e.g. V(X) = 
^ti(X 2 ). In the case of GOE, the Aijy is given by the set of N x N real symmetric matrices, and in 
the case of GSE, J\4n is the set of NxN self-dual Hermitian matrices with quaternionic entries. The 
moment matrices for both cases are given by skew-symmetric matrices and the r-functions (|3.20p 
are Pfaffians of the moment matrices [I] . In the following two subsections, we give explicit forms of 
the moment matrix and the r-functions for those models. 

3.3.1. The matrix model of GOE. The moment matrix associated to GOE is given by 

x l y°e(x - y)pt(x)p t (y) dxdy , 

where e(x) = sign(x), and pt(x) = cxp(— V(x) + 0(t, x)). 

For our finite dimensional case, fix the eigenvalues of the L of the Pfaff lattice, z± < Z2 < ■ ■ ■ < z-x n , 
set dz to be Y) S(z — Zj)dz. Then we have: 

(3.21) mtj = (zkZlT- 1 ^-* - zl-^EkEt , for 1 < i < j < 2n , 

l<fe<i<<2n 

where Ek — exp(— V(zk) + 0(t, Zk))- This is obtained from (|3.16j) with the choice of bk,i = sign(Z — k) 
for the i?-matrix (|3.11|) . The r-functions are then given by 



(3.22) T 2k = pf (M 2fc ) - > \\\\z il -z ij \]E il ---E., 



I<il<-<i2fc<2n 




J *2fc 
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One should note that the coefficient of each exponential term in (|3.22p is positive, and this gives a 
generic orbit of the Pfaff lattice on the homogeneous space SL(2n) / Sp(n) as will be discussed in 
Section 

3.3.2. The matrix model of GSE . The moment matrix associated to GSE is given by 

{z\ zi}p t (z) 2 dz = (j -i) f z*+3- 3 p t (z) 2 dz . 

I JR 

with {/, g} = fg' — gf . For the finite dimensional case, we have 

n 

(3.23) mi,j = {j 4 +j ~ 3 El , for 1 < i < j < 2n . 

k=l 

To obtain this from (|3.16[) . we first take 62fc-i,2fc = Z2k ^ Z2k \ an d — for all other (i, j). Then we 



take the limit Z2k — * Z2k-i giving doubly degenerate eigenvalues. In the formula, we have relabeled 
the indices of Zk (i.e. k for 2k — 1). Note that the matrix L in this case has n double eigenvalues. 
The T-functions are then given by 

(3.24) r 2fe = pf(M 2fe )= Y, (n^'-^n^'X- 

l<ii<— <i k <n \j<l J 

We emphasize that the form of this r-function is precisely that of the partition function of the GSE 
model. The fact that the associated eigenvalues appear as doubles is an important consequence of 
the self-dual quaternionic structure of the matrices in the GSE. 

One should compare (|3.24|) with the T-functions of the Toda lattice of st(n,R) which are given by 



Toda 




Y[\ Zll - Zl f ) E n ■ ■ ■ E lk 



We then expect that the isospectral variety of the Pfaff lattice for GSE model has a similar structure 
as that of the Toda lattice, which is that of the permutahedron generated by the orbit of the 
symmetric group S n (see for example [6j[22]). The isospectral variety for the GSE-Pfaff lattice will 
be discussed in Section [5] 



4. Real solutions 

Here we consider generic solutions of the real Pfaff lattice whose initial matrix L(0) has all real 
and distinct eigenvalues. We start with a brief summary of the results on the skew-orthogonal 
polynomials obtained in [HE] (in particular, we need Theorem 3.2 in [2]). The goal of this section 
is to determine all the fixed points of the Pfaff lattice and to find the asymptotic behavior of the 
solutions. The geometric structure of the solutions will be discussed in the next section. 

4.1. Skew-orthogonal polynomials. Let us first recall, 

LQ = QL(Q) = g* C7 7 * 1 , 

where i(0)$ = L(0)^ a V = $ A = ^IqC^V with ^ € Ge and the matrix C 7 is the companion 
matrix (|1.8|) . Recall also that the elements on the super diagonal of L are normalized to all be 1. 
The eigenvector <f> of L then defines skew-orthogonal polynomials as follows: 



L(f> = z<p, = Q* o x, 
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with x = (I, z , ■ ■ ■ , z ) ■ From </> = Q^ox with Q, ^0 € Gt, each element of 4> = {4>o> <t>i, ■ • • j <^2n-i) 
is a monic polynomial of deg(</>fc) = fc, and they are given by ipoi 2 ) — 1 ; 4>i( z ) = z an d 

2k 



(4.1) 



02fe(z) 



,2k 



l2k,2k-j 



y 2k-j 



3=1 



2k 



k = 1, . 



,n — 1 . 



<p2k+l( z ) = Z 



2k- 



1 + ^ l2k+l,2k-j z 



2fc-j 



where {(g^ 
r2 



< i , j < 2n — 1} is the matrix Q := Q'J'o. We now see that the equation QMQ T = 
H 2 J for the moment matrix gives the skew-orthogonal relation of the polynomials 4>j{ z )'- First 
define the matrix M. := V^ 1 MV~ T . Then the matrix $ = QV = ((/)i-i{ z j))i<i,j<2n satisfies 
&A4& T = H 2 J, which gives the skew-orthogonal relations, 

(<p2i,<l>2j+i)M = -(4>2j+i,<t>2i)M = h t$i,j ) < i,j < n - 1, 

and all other cases are zero. Note from M — V \MV T that the entries rriij of M give the moments 
with respect to the measure M, 



M — 



E 



i-l i— 1 



l</c,Z<2n 

where M = (mi J )i<i J <2n and .M = (/Xfc,i)i<fe,i<2n- 

Now solving the skew-orthogonal relation ^.M^ 7 = -ff 2 J, we obtain: 

Proposition 4.1. (Theorem 3.1 of PQ) XTie skew- orthogonal polynomials 4>k( z ) can be found by 







<i>2k( z ) 



T2k(t) 



Pf 



m l,2fc+l 
W2,2fc+1 



02fc+l( z ) 



72fe(t) 



Pf 



mi, 2 




m 2 ,2k 



1\ 

z 2fc 
1 



,2/c-l 



™2 : 2fc+2 

fn 2 k,2k+2 



—z 



2fe+l 



V 

for k = 1, 2, . . . , ri — 1, wra£/i <fio( z ) — 1 ^"lO 2 ) = z - 

To prove this result one notes that these functions 4>j(z) satisfy the orthogonality conditions, 
(fok, zJ )m = , for < j < 2k, and (0 2 fc+i , z j ) M = , for < j < 2k - 1 and j = 2k + 1. 
The following forms of the skew-orthogonal polynomials are particularly useful: 

Proposition 4.2. (Theorem 3.2 in [2]) The skew- orthogonal polynomials 4>k{ z ) can be expressed in 
terms of r- functions, 

'<hk(z) = -^T M (t- [z- 1 ])**, 

1~2k\t) 



[0 2fc+1 (z) = -i^(z + A) 7afc (t-[z^])z» 
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where r 2k {t - [z 1 ]) = T 2k {ti -\M~ ■■■)■ 

(A proof of this Proposition is based on the moment equation given in (|3.17|) . and Lemma l4~T1 below.) 

From Proposition 14. 2\ we have the explicit form of the entries qij of Q in (|4.1|) in terms of the 
r-functions, 



(4.2) 



Q,2k,2k-j{t) 



s 3 (-d)T 2k (t) 



Q2k+l,2k-j(t) — 



T2k (t) 

[s j+1 (-d) + d 1 s 1 (-d)}T 2k (t) 

T2k(t) 



Here Sk(—d) denotes Sk(— £?i, — \di, — §03, • • •)) w ith dk = and Sk(ti,t 2 , . . .) are symmetric 
functions defined by 

(oo \ oc 

j=l J fc=0 

The explicit form of the Sfc(ii, . . . , t&) is given by 



Sk{ti,t 2 , . . . ,t k ) = \ ~~T~~y —. , 



'_k_ 

ji+2j 2 +-+kj k =k Jllj2] ' ' ' Jk ' 



which is the complete homogeneous symmetric function hk(x\, . . . , x 2n ~i) with ife = r 53i=i 1 a;fe 
(see [20]), 



«5fe(tl, . . . , t k ) = hk(xi, . . . , X2n-l) = 



X 2 J-2n-l 



il+«2H M2n-l=fc 

With the Q given in (|4.2p . we can find L in terms of the r-functions, 

(4.4) L = QC^Q- 1 . 

In particular, the 2x2 block diagonal part of L has the form, Lq := diag 2 (£o,o, £i,i> . . . , -L rl _i. n _i) 
with 

(4.5) L k , k =( q2k ' 2k - 1 1 ) fc = l,2,...,n-2, 

\92fe+l,2fc-l — <72fc+2,2/c ~92fe+2,2fc+l / 



and 



7 ■ _ / 1 \ i / <&n-2,2n-3 1 

( ^92,0 -92,1/ V 92n-l,2n-3 ~ 72 



Notice from (|4.2p that the diagonal elements of L (also for L in (| 1 . 7[) ) are given by bk — — Q2k,2k-i = 
hi T2fe . The formulae for bk together with a k in (|3.8p are similar to the formulae for the dependent 
variables in the Toda lattice, i.e. 

, . ^T 2k - 2 T 2k+2 d 

(4.6) a k = , bk = ^—lnT 2k , tor k = 1, ...,n-l. 

T 2 k oti 

In particular, the functions bk(t) determine the diagonal elements of L, and provide the basic infor- 
mation about the geometry of the isospectral variety, similar to the Schur-Horn theorem for the set 
of symmetric matrices (see Section [5]). 
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4.2. Fixed points of the Pfaff flows. We now compute the asymptotic values of L for the case 
of L(0) whose eigenvalues are all real and to be ordered as 

Zl < z 2 < ■ ■ ■ < z 2n . 

In the case of GSE-Pfaff lattice, we take the limit z 2k — ► z 2 k-i- 
Let us start with the following Lemma for computing Sj(—d)T 2 k: 

Lemma 4.1. With Ek{t) = exp9(t, z k ), we have, for m > j , 

Sj{-d) E 1 --E m = (-l)*aj(zi, ...,z m )Ef-E m , 

where o~j(z\ . . . , z rn ) is the elementary symmetric polynomial of degree j in (zi, ... , z m ), 

^ j (zi j • • • i Zrn) — ^ %i\ ' ' ' Zij . 

l<ii < ■ --<ij <m 

(note 7fe = crfe(zi, . . . , z 2n ) in C 1 ). For j > m, we have Sj(~d)Ei ■ ■ ■ E m = 0. 
Proof. Since -§t~E k — J^r-E'fc, we have 



^-[^])=exp[-|:^A 

exp In 1 — 

V z oti 



, ^ 1 03 

EJt) = exp - > — r — , 
1 f^Jzidt{ 

Eiil) = ( I -^)E k (t). 



E k {t) 



Then from f|4. 3|) we find 

-S 3 (-d) {Eh.-. E m )(t) = (E 1 -. E m )(t - [z- 1 ]) 



3=0 



n(i-^) (E 1 ...E m )(t). 



i=l 



Comparing the coefficients of z~ k , we obtain the desired formulae. □ 

This lemma confirms again that 4> 2n {z) = det(L — zl). From this Lemma and the ordering of the 
Zfe's, we have, for t\ — > oo, 



(4.7) 



q2K2k- ] {t)^{-iycj ] {2k), 

«a*+i,2*-j(t) — > (-iy +1 a J+1 (2k) + (-iya 1 (2k)a j (2k) . 



Here we introduced the notation cr 1 (2fc) :— cr 1 (z 1 , . . . , z 2k ). We then have 

Theorem 4.1. In generic case, the L(t) approaches the block upper triangular matrix L~ as t\ 



( L 0,0 f" °2 

2 e_ 



L~ = 



o 2 o 2 ••• l: 



V o 2 o 2 



o 2 \ 
o 2 



e So © 0+ , 
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where e_ = ^ jjj and X fc fe is the 2x2 matrix defined by 

L- -( - ai(2k) 1 W-0 1 n-1 

k ' k ~ {- (o-!(2k) + z 2k+1 ) (a^k) + z 2k+2 ) a 1 (2k + 2)J A, 

wit/i <Ti(m) = (7i(zi, . . . , 2 m ) = X3j = i 2j and <7i(2n) = tr(L) = 0. 

On i/ie other hand, as t± — > oo, the L(t) approaches L + whose 2x2 block diagonal elements are 
given by L+ = ,_ 1 for k = 0, 1, . . . , n - 1. 



Proof. First note from (|4.2p that Q = Q^o approaches a constant matrix as t\ — ► ±oo. This 
implies Q _1 f^ — » as t\ — » ±oo. On the other hand, from Lemma f3. H we have 

B X = H (L- - JL\.J) H~ l =L-- JHI^H^J . 

The non zero entries in H(L + ) T H^ 1 are given by a\ (see (|3.8p ). With the ordering of the eigenvalues, 
i.e. z\ < z 2 < ■ ■ ■ < z 2n , and the genericity, one can see 



T 2k {t) 



\c Ei---E 2k , as ti — > -oo; 

\c + E 2n -2k+i ■■■E 2n , as t\ > oo, 

where are nonzero constants. Then from ()3.8j) . we have a k (t) — > as t\ —> ±oo, therefor L + — ► 0. 
Since Bi — > and Jf/X+i? -1 J — > we conclude that Z_ — > 0. 

Since L + is just a constant matrix, it only remains to find the limits of Lq as t — > ±oo. The Lq 
is explicitly given by (|4.5[) . and using (|4.7[) . we can easily compute Lq for t\ — * ±oo. □ 

Each submatrix LjT fc in the asymptotic value of L has the eigenvalues {z 2k +i, z 2k+2 }. This 
implies that the Pfaff lattice has the sorting property of eigenvalues in pairs as t\ — > ±oo, i.e. 
({zi, z 2 }, . . . , {z 2n -i, Z2n}) as ii — * -oo and ({z 2n ^i, z 2n }, . . . , {z iy z 2 }) as <i — > oo. It should be 
noted that the order in the pair is free. 

We can also show: 

Lemma 4.2. 4 matrix L S Zj&(7) is a /tied point of the Pfaff lattice, if and only if L(t) is a block 
upper triangular matrix for all t. 

Proof. Let us first note that the Pfaff lattice equation for L is given by 

dL 

where Bj is given in Lemma 13.11 In particular, we have B\ = L_ — JHL^H^J, and with a k = 
hk+i/hk, the second term can be explicitly expressed as 

n-i 

JHL 1 j r H~ l J = JH 2 L+H~ 2 J = ^2 a k e 2(fc+l),2/c-l ; 

fc=i 

where e^j is the 2n x 2n matrix with 1 at (i, j)-entry and for all others. Then if L is block upper 
triangular, i.e. £_ = 0, we can compute directly that is strictly block lower triangular with 

entries given only by the a^'s. This implies = 0, and therefore = which fixes a k = 

for all k. One also notes that if L(t) is not block upper triangular, then ^ 0, therefore the only 
possible fixed points are those of block upper triangular shape. A similar direct calculation gives 
the same conclusions for the higher flows with generator Bj for j > 2. □ 

We can also obtain the explicit block upper triangular form of a fixed point, L £ Z^(j), of the 
Pfaff flow; we find that a general fixed point is in a form similar to that of L~ in Theorem 14. II 
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t2 - t2 

50 




-80 -40 

t3=10 



40 t 1 



Figure 2. The evolution of &i(ii, ti, t%), the diagonal element of L for the GOE- 
Pfaff lattice with n = 2. The pair of numbers {i,j} indicates the value of b\ in the 



corresponding region, i.e. b\ 



j. Notice that the pairs {1, 3} and {2, 4} appear 



Lk,k 



i the interaction region forming a double quadrangle cone in the ^-direction with 
a singular point at t\ = t 2 = ts = 0. 

Theorem 4.2. Each fixed point of the Pfaff lattice is uniquely parametrized by a set 

({z h , z l2 },..., {z i2n _ 1 , z l2n }) : with i 2k -i < i2k, k = l,...,n. 

Therefore the total number of fixed points is given by (2n)!/2™. Moreover the fixed point corresponding 
to this set is expressed by a block upper triangular matrix in Z^{^) whose block diagonal is given by 

diag 2 (i) = (i ,o,£i,i, • • • ,L n -i,n-\) with 

-ai(ii,...,i 2 k) 1 

-(<7l(n, . . .,i 2 k) + Z l2fc + 1 )((Tl(n, • • .,12k) + Zi 2k+2 ) <7l(il, . . .,i 2 k+2) / 

where t7i(ii, ...,i m ) = J^jLi z h and CT i( ) = °- 

Proof. For a block upper triangular matrix with the eigenvalues {zi, . . . , Z2 n }, each 2x2 matrix 
in the 2x2 block diagonal part has a pair of eigenvalues {zt 2k _ 1 , Zi 2k } with 1 < i 2 fc— l < *2fe < 2n, 
therefore a fixed point of the Pfaff lattice gives the set above. Then from the form of L in Z^j), one 
can construct a unique L corresponding to this set: In the first block (the top left 2x2 submatrix 
Lo,o), w e have a unique block, 

1 

-Z\ 1 Z{ 2 Zj jl -\- Zi 2 ^ 

Now suppose L>k, k has the above form, then, as L £ ^(7), ifc+i,fc+i has the form 

-0-\{h, . . . ,l2k+2) 1 

x b k+2/ 

Since this block has the eigenvalues {zi 2k+3 , Zi 2fc+4 }, one finds the unique bk+ 2 — —o-i{h: ■ ■ ■ , hk+4) 
and x = -(<7i(ii, . . . ,i 2k +2) + 2i 2H3 )h(ii, • • - ,^+2) + Zi 2k +J- Tm s completes the proof. □ 

In the next section, we identify those fixed points as the vertices of the moment polytope generated 
by the flows associated with the Chevalley invariants, that is, the torus-fixed points of the Pfaff 
lattice. 

Example 4.3. The case n = 2 for the GOE-Pfaff lattice: The T2-function is given by 

l<i<j<4 



GEOMETRY OF THE PFAFF LATTICES 



2!) 



where E(i,j) = (zj — Zi)E{Ej with the ordering z\ < z 2 < z 3 < 24. We have 4!/2 2 = 6 fixed 
points, and each fixed point corresponds to the exponential term bijE(i,j) for 1 < i < j ' < 4. 
We consider the Pfaff orbit in the space of (ti,t 2 ,t 3 ) corresponding to the Chevalley invariants, 
Hk(L) ■= jrpjtr(X fe+1 ) for k = 1,2,3, that is, we consider (|1.6|) . We consider the asymptotic values 
of the diagonal element of L(t), i.e. 

6i(*) = ^-ln75(t). 

If the exponential term E(i,j) in the t 2 dominates the others in some region of (ti,t 2 ,t 3 ), then in 
this region b\ takes the value 

&l(tl)*2) £3) ~ z i + z j ■ 

One can then identify each fixed point with the ordered set of pairs [{zi, Zj}, {zk, zi}) : for i < 
j, and fc < Z as the region where E(i,j) gives the dominant term. So we expect to have six different 
regions in {ti,t 2 ,t 3 ) space: Figure [2] plots b\ (t) for the case of n = 2 in the GOE setting where 

fO 1 1 1\ 
Oil 

B - 1' 

V 0/ 

with (zi, Z2, z 3 , Z4) = (—2, —1, 0, 3), that is, we have 

r 2 (t) = E1E2 + 2E 1 E 3 + 5£i£ 4 + E 2 E 3 + AE 2 E A + 3E 3 E± . 

In Figure [21 each region with the dominant E(i,j) is marked by {i,j}, and the boundaries of 
the regions are given by the balance of two exponential terms, e.g. E\E 2 = 2E\E 3 which gives 
9(t, z 2 ) = 9(t, Z3) + In 2. This equation describes the boundary line between {1, 2} and {1, 3} in the 
ti-t 2 plane for each fixed t 3 , 

(z 2 — zi)ti + (z| — z\)t 2 + (z| — zf)t 3 — constant. 

This line expresses a line soliton solution of the coupled KP (or DKP) equation, and Figure 
illustrates a resonant interaction of 2-soliton solution (see |19|). 

4.3. The GSE-Pfaff lattice. Recall that the moment matrix for the GSE model was obtained by 

taking the limits z 2 k — > z 2 k-i, and setting &2fc-i.2fc = — r 24 for some constants Ck and bij = 

for all others. Note here that if Ck > 0, then it can be absorbed into the exponential term E 2 k-iE 2 k- 
However those constants will be important for the general solutions including rational solutions of 
the Pfaff lattice (see below). We also relabel the eigenvalues z 2 k~i as Zk for convenience as in Section 
13X21 Then we have 

Theorem 4.4. As t\ — > —00, the L(t) for the GSE-Pfaff lattice approaches the upper triangular 
matrix L~ , as in Theorem \4.1\ whose block diagonal elements are now given by 

f- _( -2(Ti (fc) 1 

- y-(2at(k) + z k+l f 2a x (k + 1) 

with CTi(fc) = z j> an d a i(®) = cr i( n ) = 0- 

Proof. In the limit z 2k — > ^2fe-i, the o"i(2fc) becomes 2 S^=i z ^j-li which is just 2<7i(fc) after 
relabeling the eigenvalues. □ 

Notice that the matrix L^ k has doubly degenerate eigenvalues Zfe+i,. Proposition 14.41 means that 
the GSE-Pfaff lattice flow sorts the double eigenvalues in the same manner that the Toda lattice 



0,1, 
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Figure 3. The graphs of the diagonal elements b\ and b 2 of L for the GSE-Pfaff 
lattice with n = 3. The number sets {i} and {j, k} with j < k indicate the values 
of bi and b 2 , i-e. b\ — 2zt and 62 = 2(zj + Zk) in those regions. 



flow sorts eigenvalues. One can also see that there are n\ number of fixed points, and each fixed 
point corresponds to the matrix having the permuted double eigenvalues, 

( z ii ? z i 2 ? ■ • ■ : z i„ ) ■ 

Those points form the vertices of the permutahedron of the orbit of the permutation group, which 
is the isospectral manifold of the Toda lattice (see e.g. [22j[6], also see Section l573|) . 

Example 4.5. The GSE-Pfaff lattice with n = 3: The r-functions are given by 
r 2 = Cl E\ + c 2 E\ + c 3 El , 

r 4 = c lC2 (zi - z 2 fE\E\ + cxcsiz! - z 3 fE\E\ + c 2 c 3 {z 2 - z^E\eI , 

where are nonnegative constants, and E^ = e e ^ t,Zk ^ with 6{t,Zk) — Y^j=i^j z l- Each exponential 
term in r 2 k becomes dominant in a certain region of the (ti, t 2l £3) space. Figure [3] plots the functions 
bk = ^-lnr2fc for k = 1,2 to show the dominant exponentials in the r-functions. Then we obtain 
the fixed points of the Pfaff lattice. For example, the region having b\ = 2z 2 = {2} and b 2 = 
2(zi + z 2 ) = {1, 2} gives the fixed point L~ = diag 2 (Lo~ , X -E 22 ) with 

\\ f _ f -2z 2 1 \ ~ {2z 3 1 N 



^o,o - 



4 2z 2 

where we have used z\ + z 2 
(z 2 ,z 1 ,z 3 ). 



-(Z2-Z3) -2z 3 







Z3 



0. This fixed point corresponds to the ordered eigenvalues 



Remark 4.6. The rational solutions of the Pfaff lattice discussed in [3] are given by the GOE r- 
functions in the nilpotent limit with specialized choices of the i?-matrix. (By the nilpotent limit, we 
mean the limit Zk — > for all fc.): We compute the moment matrix M — EBE T = e e ^' C ^ Be 6 ^' ^ 
with B = VBV T in the nilpotent limit corresponding to 7 = (i.e. Zk = for all fc). Choosing an ap- 
propriate B-matrix so that there exists a limit B$ — lim 7 ^o VBV T , we have Mq = e 6 ^' 00 * 1 Boe 8 ^' ^ 
with 

(1 5i(t) S 2 (t) ••• S 2 n-l(*)\ 

1 ••• S 2n -2{t) 

J(t,C ) = 1 ••• S 2n -3(t) 



\0 



1 / 
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where Sk(t) = Sk(ti, . . . , tk) are the complete homogeneous symmetric functions denned in (|4.3I) . 

Various choices of Bq are considered in [3] , some of which contain the Jack polynomials. In particular, 

if i?o is chosen as the matrix having only nonzero anti-diagonal elements with &fc,2n-fe+i = 2(n — 

fc) + 1 for fc = 1, . . . , n, the Pfaffian T2k — pf((Mo)2fc) then gives the Jack polynomial J^ 2 \x) for 

fc 

rectangular partitions, A = (2(n — fc), . . . , 2(n — k) ) with tj = i a^', which is a zonal spherical 
function on GL(2n)/Sp(n). It is then interesting to compare this with the case of the nilpotent Toda 
lattice of sl(n) where the r- functions are given by the Schur polynomials with rectangular partitions, 

fc 

i.e. T l oda = (-l) M V 11 5A(ti,. .. ,i„_i) with A = (n — fc, . . . , n — fc) for fc = 1, . . . , n — 1 (see [7]). 

5. Geometry of the isospectral variety 

Here we discuss the geometric structure of the isospectral variety for the Pfaff lattice using the 
moment map defined in [10) . 

5.1. Geometric structure of the r- functions. Recall that the solution of the Pfaff lattice is 
given by the matrix factorization, e e( *' i(0)) = Q{ty 1 P(t) and L(t) = Q(t)L(0)Q(i) _1 . This implies 
that we have a bijection between the isospectral variety Zr(7) in (|1.9[) and the homogeneous space, 

SL(2n,R) 
■ Sp(n,R) ' 

with dim7i = 2n 2 — n — 1. Note that H is a symmetric space associated with the involution 
6»(expX) = exp(<r(X)) on SL(2n) with a(X) = JX T J, i.e. a = d9. Then each point of the Pfaff 
orbit can be parametrized by a skew-symmetric matrix. The B-matrix used in the r-functions indeed 
gives a parametrization of the point of H.. Namely we have 

Zr{j) — r -+ Tt Skew{2n) 

L(0) i — > ^o" 1 modSp(n,R) i — > B = ^ 1 J^- T 

where B = VBV T and L(0)^q = ^oC-y with £ Gj. Here the <j) gives an isomorphism when <fi is 
restricted on an open dense subset of Tl, and Skew(2n) is the set of skew-symmetric matrices with 
dct = 1. The _B-matrix then parametrizes an initial point of the Pfaff flow, and the time evolution 
follows 

L(t) i — ► e «(*>c T )^-i mo d5p(n,M) i — ► M = E(t, A)BE(t, A) T 

Note here that M = e e(t ^Be e ^ c ^ T with E(t,A) = Ve 6 ^ = e e ^ c ^V . The Pfaff flow also 
defines a torus action on TL, that is, e 6 ^' ^^ 1 Sp(n) = Ve 6 ^^^ 1 Sp{n). 

Each r-function T2fc can be expressed in the following form, which is a Pfaffian version of the 
Binet-Cauchy theorem |14) : 

Lemma 5.1. 

T 2 fc = pf(M2fc)= P f (*i: ■ ■ • ,hk)E(ii, . . . ,i 2 k) ■ 

l<ii<--<t2k<2n 



where 



:S2 
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• E(ii, . . . , i2k) is the determinant of the 2k x 2k submatrix of the Wronskian matrix E(t, A) 
in 13.12\) . and is given by 



E(h,. 


■ ,hk) = det 


K 


Ei 2 

K ■ ■ 


■ E i2k \ 

■ E L 






Uf- i} 


E i2 


P (2fc-1) , 



• pf (ii, . . . , i 2 k) is the Pfaffian of the 2k x 2k skewsymmetric submatrix of the B-matrix, and 
is defined by 



pf(«l, . . . ,12k) = Pf 



/0 b iui2 




V 



l 2k-l-,' l 2k 





Recall that {Ei : i = 1, . . . , 2n} forms a basis of K 2 ™ for distinct Zk, let us further assume that 
z h,...,i 2 k — z ii + ' ' ' + z i-2k are au distinct. There is a bijection between E(ii, . . . ,i 2 k) and the basis 
vector ei 1 A ei 2 A • • • A ei 2k of A 2fe M 2n , but this bijection is not a homomorphism of the A-product 
structure on the space ©jLj A 2fc M 2 ™, this deficiency does not affect our computation. The coefficient 
of the E(ii, . . . , Z2fe) gives a Pfaffian-Pliicker coordinate, which is written recursively by the expansion 
formula, 

2k 

pf(ii,i 2 , . . .,i%k) = E' : 1 : '''< • P f (* 2 ' ••■>h>-- ->*2jfe) > 

where ij implies the deletion of ij from pf(ii, . . . , «2fc), and pf(z, j) — bij. For example, the Pfaffian 
pf(l,2,3,4) is given by 

pf(l, 2, 3, 4) = 6i, 2 6 3 ,4 - &1,3&2,4 + &1,4&2,3 • 

One should also note from (14. 2[) that each element q r ^ in Q is invariant under the scalar multiplication 
to T2fc as is L from (|4.4[) . that is, T2fc = c/.T2fe for any constant Cfe. This means that the T2fe can be 
considered as a point of the projective space P(A 2fe Kl 2 ™). 

We note here that T 2 k is not (in general) decomposable in A 2/c K 2n , therefore T 2 k does not represent 
a point in the Grassmannian Gr(2k, 2n). 

5.2. Moment polytope. Let us recall that the unipotent subgroup G{ can be identified as the top 
cell Y^id of the flag variety SL(2n, M.)/P 2 x2 (see Section [3~TT) . It is also well-known that there is an 
embedding of the flag SL(2n,R)/ P 2x2 , 



SL(2n,R)/P 2x2 ^ P(A 2 M 2 ") x P(A 4 M 2 ") x • • • x P(A 



2n — 2Trj)2n 



This suggests that G{ can be embedded in this product space (see (I3.6[) V 

On each P(A 2fe R 2 ™), one can define a moment map, \i : P(A 2fe M 2 ™) — ► f)Jj, with a torus action 
given by the Pfaff flow g — e e (*' A ) = diag(£i, . . . , E 2n ) (see [TUJ [3T]): By considering r 2k to be a 
point on P(A 2fe M 2l; 



(5.1) 



), we define the image fJ,(r 2 k) as follows: For the expression of r 2 k in Lemma 15.1 



H{ T 2k) = 



E ip f ^< 

\<ii<---<i 2k <2n 



,i2k)E tl ■ ■ -E i2h \ (Ci ± 



E 

l<il<-<i 2 k<2n 



|pf(ii 



,i 2fe )^---P i2 J 2 
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Here f)Jj is the dual space of Cartan subalgebra of st(2n, R), 



()r = Span K < £1, . . . ,C 2n : ^ A 



j>2n— 1 



with the weights Ck- The moment map [i can be extended to the map over the product space of 



» 2fciTj2r 



)by 



(5.2) 



2n ) x 
(t 2) 



x F(A 

! T 2n-2 



2n — 2TiT)2n\ 



M T 2fc) • 



The moment polytope given by the image of the moment map is expressed by the wedge product 
representation with the highest weight (see e.g. [TT]). 

n—l n — 1 

w* ^ (Ci H h £ 2 fc) = X]^ n - fc )( £ 2fe-i + ^2fe) • 

k=l k=l 

This is a tensor product of the basic representations of SL(2n), and it is customary denoted by 
To, i,o, ...,i,o ( see [H]) : In general, r aii .... a2 „_i represents the irreducible representation of SL(2n) 
with the highest weight w* = X^fcHi 1 a k{£i + ••• + •£*). The vertices of the polytope are given by 
the weights which are parametrized by the set of numbers {0, 1, . . . , n — 1}, i.e. 

2/i 

w a = (a%, a 2 , . . . , a 2 „) := ^ ak ^ k • 

fe=i 

with ckfc € {0, 1, . . . , n — 1}, and each number of the set appears exactly twice in e.g. if* = 

(n— 1, n— 1, n— 2, n— 2, . . ., 1, 1, 0, 0). Each vertex of the polytope can be parametrized by a unique 
element of W p , the set of minimal representatives of the cosets in, 

W p := S 2n /Wp with W P = (si, s 3 , . . . , s 2 n-i) , 

where Sfc are the simple reflections, i.e. Sk — (k,k + 1) for k = 1, . . . , 2n — 1. Note here that 



Wp = S 2 x • • • x S 2 and |W P | = (2rt)!/2™. For example, for n — 2 and Wp = (si, S3), we have 

W p = { e, s 2 , sis 2 , s 3 s 2 , S1S3S2, S2S1S3S2 } • 

We also note that the quotient W p parametrizes the Bruhat decomposition of SL(2n)/ P 2x2 , i.e. 
SL(2n)/ P 2x2 = U weW p NwP 2x2 where N is a unipotent subgroup of lower triangular matrices. 

From (|5.ip . one can also see that those vertices are identified with the fixed points of the Pfaff 
lattice. The identification is given as follows: For the vertex of (a\, . . . , a 2n ), assign o^'s so that 

n — fc— 1 k = 0, 1, . . . , n — 1 . 



a 



L t2k + 2 



Then the corresponding fixed point is given by the block upper triangular matrix L obtained in 
Theorem 14.21 For example, the vertex w a = (0, 1, 2, 1, 0, 2) in the case of n = 3 corresponds to the 
matrix L given by 

/ 1 \ 

-z 3 z 6 z 3 + z 6 1 

-(z 3 + z 6 + z 2 )(z 3 + z 6 + z 4 ) -(zi+z 5 ) 1 

z\ + z 5 1 

\ -zxz 5 0/ 



L = 



34 



YUJI KODAMA* AND VIRGIL U. PIERCE" 




(0110) 1 1 

Figure 4. The momentum polytope for the GOE-Pfaff lattice with n = 2. The 
right figure shows the function b\ of the Pfaff lattice in (ti,t2,ts) with i 3 = —10 
(see Figure [2]) . Each asymptotic region of the b\ is identified by the vertex of the 
polytope. 

where we have used tr(L) = ai(l, . . . ,6) = 0, and all other entries are zero. This fixed point can be 
parametrized by the element S1S2S1S3S5S4S3S2 € Se/(s±, S3, S5). 
We summarize the result in this section as follows: 

Theorem 5.1. For the generic Pfaff lattice, the image of the moment map i f 5. £j) is a convex polytope 
whose vertices are parametrized by the set W p of minimal representatives of the cosets in Si n jWp 
with Wp = (si, S3, . . . , S2n-i). Those vertices are identified as the fixed points of the Pfaff flow given 
in Theorem \4-2\ 

Example 5.2. In Example 14. 3 [ we considered the Pfaff lattice for the GOE model in the case of 
n = 2. The momentum polytope associated to this model is given by the irreducible representation 
with the highest weight w* = C\ + C2, which is an octahedron. In Figure^ we show the momentum 
polytope and explain how the polytope can be realized from the solution of the Pfaff lattice in the 
space of (t\,t2,tz). One should note that the polytope is just a dual diagram of the pattern of the 
solution b\ in (£1,^2,^3), that is, identify each region for constant b\ as a vertex of the polytope, 
and the line between two values of b\ as an edge of the polytope, and so on. The vertices are 
parametrized by the elements of W p , i.e. 

e={l,2}, s 2 = {l,3}, Sl s 2 = {2,3}, s 3 s 2 = {1,4}, 
S1S3S2 = {2,4}, S2S1S3S2 = {3,4}. 

The generic orbit can be described by a curve inside of the polytope approaching to {1,2} vertex 
as ti — * —00 and to {3, 4} vertex as t\ — > 00. Those vertices represent the highest and lowest weights. 
The edges of the polytope correspond to non-generic orbits of the Pfaff lattice. For example, the 
edge between {1, 2} and {1, 3} is a solution given by the r 2 function with the B-matrix, 

/0 1 1 0\ 


B - oo- 
V 0/ 

Note that the T4 = for this choice of _B-matrix, and the T2 satisfies the KP equation. The faces 
of the polytope correspond to non-generic orbits (i.e. KP solutions). Some explicit solutions are 
studied in pjj] . We then note that a classification of the solutions of the Pfaff lattice can be obtained 
by the representations of the SL(2n) orbit in P(ro 1 i,o,...,i,o)- 
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5.3. The moment polytope for the GSE model . Recall that the r-functions for the GSE-Pfaff 
lattice are given by p.24j) . i.e. 



where A(ij 



with 



l<ii <---<ifc<n 

= Y[i<j<i<k( z h ~ z ij)- Writing those in the form, 
J2 A(i 1 ,...,i k ) 3 F(i 1 , 

l<ii<-<ik<n 



■2 

ik ' 



l<j<l 
T2k 



,lk) 



F(i u ...,i k ) :=2- M ^det 





F 


F' 


F' 


n 




(fc-i) 


F (k-1) 



I 



where Ff. = E\ = e 2e ^' Zfc - ) and — d^F^fdt\. Since F(i±, . . . , ik) represents a coordinate for 
A fc R", the T2fe can be considered as a point in P(A fc M") as discussed in the previous section. Now 
we consider the moment map \x : P(A fe R n ) -> f)Jj where JjJ is the dual space of Cartan subalgebra of 

\A( iu ..., lk fF 1 ---F k \ 2 (C ll +--- + C lk ) 

I \ \<%\<—<ik<n 
M T 2fcJ = ^ " " > 

£ |A^,...,z fc ) 3 iV^.| 2 

where Ck are the weights of f)* of sl(n,M). As in the previous case, the moment map \jl can be 
extended to the product space, 



(5.3) 



x P(A 2 R") x ... x P(A"- 1 R") 
(r 2 , r 4 , . . . , r 2n _ 2 ) 



fc=l 



Then the moment polytope is the permutahedron of S n (denoted as Ii,. i of the basic representation 
of SL(n) [H]), which is the weight polytope with the highest weight, 



£(A + ••• + £*) = (n-l,n-2,...,l,0) 



fe=i 



Each vertex of the polytope is parametrized by the element of S n . We summarize the result: 

Theorem 5.3. For the GSE-Pfaff lattice, the image of the moment map H5.3\) is given by the 
permutahedron of the symmetric group S n , whose vertices are the fixed points of the flow. 



Example 5.4. The momentum polytope for n = 3 GSE model: In Figure El the dominant expo- 
nential terms in the r-functions are plotted in plane. The combined graph of b\ — <91iit 2 /<9<i 
and £> 2 = dlnr^/dti divides the £i-t 2 plane into 6 regions, each of which is marked by the number 
set ({i},{j,k}) with i = j or % = k (j < k) (see Figure [5|). The dual graph of the combined graph 
gives a hexagon as shown in Figure [5j which is the moment polytope for the symmetric group S3. 
The edges of the hexagon correspond to non-generic flows, e.g. the edge between the vertices (201) 
and (102) corresponds to the flow given by the r-function with c 2 = 0, 

T2 = ctEl + c 3 El . 
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Figure 5. The momentum polytope for the GSE-Pfaff lattice with n = 3. The right 
figure shows the combined graph of bi and 6 2 in Figure[3] Each ({i}, {j, k}) indicates 
the set of asymptotic values (61,62), i.e 61 = Zi and 62 = Zj + Zk- The left figure 
shows the moment polytope with the weights (a\, a%, a.3) — aiCi+a2C2+<^3J~-3, and 
«fe indicates the number of k appearing in the set ({?}, {j, I}). The Weyl elements 
Si,S2 and S1S2S1 correspond to the boundaries of the regions ({«}, {j, k}). 

This then gives r 4 = 0103(2:1 — z 3 ) 4 EfE^, and the Pfaff lattice is reduced to the model with n = 2, 
i.e. a subsystem of the original lattice. 
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